LIBRARY 
«VAt  POSTGRADUATE 
MONTEREY,  CA! 


ft^aifffnkL 


t.fw^' 


IS 


CONFroENTIAL 


A  THEORETICAL  STUDY  OF  AUTOMATIC  INERTIAL  NAVIGATION 


Leonard  Erb 
Lewis  J.  Stecher,  Jr, 


by 

BJ3.,  United  States  Naval  Academy,  1941 
BJS.,  United  States  Naval  Academy,  1941 


SUBMITTED  IN  PARTIAL  FULFILLMENT  OF  THE 

REQUIREMENTS  FOR  THE  DEGREE  OF 

MASTER  OF  SCIENCE 

at 

MASSACHUSETTS  INSTITUTE  OF  TECHNOLOGY 

1949 


■■ 


• 


UMMeMM 


This  document  contains  information  affecting  the  national  defense  of 
the  Unitea  States  within  the  meaning  of  the  Espionage  Act  50  U.S.C., 
31  and  32  as  amended.  Its  transmission  or  the  revelation  of  its  con- 
tents in  any  manner  to  an  unauthorized  person  is  prohibited  by  law. 


™*fi&tm\Ai 


May  20,  1949 


Prof.  Joseph  S.  Newell 
Secretary  of  the  Faculty 
Massachusetts  Institute  of  Technology 
Cambridge  39,  Massachusetts 

Dear  Professor  Newell: 

In  accordance  with  the  regulations  of  the  faculty, 
we  hereby  submit  a  thesis  entitled,  A  THEORETICAL 
STUDY  OF  AUTOMATIC  INERTIAL  NAVIGATION  in  par- 
tial fulfillment  of  the  requirements  for  the  degree  of 
Master  of  Science. 


tgijteslFffiD 


ACKNOWLEDGMENT 


The  authors  wish  to  express  their  deep  appreciation 
to  Dr.  C.  S.  Draper  and  to  Dr.  W.  Wrigley  for  the  interest, 
confidence,  and  assistance  given  by  them  in  supervising  the 
progress  of  this  work. 

Invaluable  aid  was  also  given  by  the  staff  of  the 
Instrumentation  Laboratory,  particularly  of  the  "Febe" 
Section.    Individual  mention  is  foregone  only  because  of  the 
many  members  who  helped. 

The  staff  of  the  Rockefeller  Analyzer,  at  the 
Massachusetts  Institute  of  Technology,  by  their  interest  and 
cooperation,  contributed  greatly  to  the  work. 

Thanks  are  given  to  Mr.  L.  E.  Payne  and  his  associ- 
ates of  Jackson  &  Mo  re  land  for  their  part  in  the  preparation 
of  the  text  and  illustrations. 


iv 


awossffiSR 


ABSTRACT 

The  purpose  of  this  thesis  was  to  postulate  a  method  of  long-range 
inertial  navigation,  using  only  acceleration  inputs,  and  to  examine  theo- 
retically the  response  of  this  system  to  external  disturbances,  using  vari- 
ous types  of  mechanizations.    The  problem  was  divided  into  one  of  track 
control  and  of  range  indication,  using  a  great-circular  path  between  the 
points  of  departure  and  destination. 

The  response  of  the  track  control  system  was  examined  with  mecha- 
nization equations  of  as  high  an  order  as  the  fifth,  and  was  found  to  improve 
as  the  order  of  the  mechanization  equation  increased.   Response  to  impulse 
type  wind  acceleration  disturbances  was  satisfactory,  but  long-period  sinu- 
soidally  varying  winds  caused  excessive  errors.    It  appears  probable  that 
additional  feedback  loops,  in  conjunction  with  a  mechanization  equation  of 
the  fifth  or  higher  order,  will  solve  the  track  control  problem  satisfac- 
torily. 

In  the  range  indication  system,  it  was  found  to  be  impossible  to  re- 
move any  forcing  function  terms  of  higher  order  than  the  time  rate  of 
change  of  acceleration.    It  appears  probable  that  the  range  indication  sys- 
tem, as  postulated,  will  prove  to  be  satisfactory,  especially  if  the  sensi- 
tivities are  made  variable,  as  a  function  of  the  angle  of  the  input  pendulum. 
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CHAPTER  I 
INTRODUCTION 

The  importance  of  developing  a  long-range  self-contained  system  of 
automatic,  mid-course,  navigation  for  the  control  of  guided  missiles  can- 
not be  overestimated  (Fig  1*1).   As  long  range  flights  of  piloted  aircraft 
over  enemy  territory  become  increasingly  hazardous,  because  of  the  neces- 
sity of  travelling  without  fighter  protection  toward  heavily  defended  cities 
provided  with  elaborate  warning  nets,  automatic  navigation  can  provide  a 
method  of  delivering  bombs  without  risking  the  lives  of  highly  trained 
pilots. 

With  the  development  of  long-range  jet-propelled  missiles,  automatic 
navigation  should  eventually  permit  the  United  States  to  bomb  any  place  in 
the  world  from  bases  located  within  its  territorial  limits. 

For  such  a  system  of  automatic  navigation  to  be  successful,  it  must 
have  a  degree  of  accuracy  comparable  to  that  achieved  by  human  pilots.   It 
must  not  be  easily  susceptible  to  "jamming"  by  enemy  action.    It  must  op- 
erate successfully  at  supersonic  velocities,  and  in  the  upper  reaches  of  the 
earth's  atmosphere. 

At  present,  the  problems  of  automatic  long-range  navigation  are  under 
consideration  by  several  groups.   Among  these  are: 

1.  Baird  Associates,  Inc. 

2.  Hughes  Aircraft  Co. 

3.  Instrumentation  Laboratory,  M.I.T. 

4.  Kollsman  Instrument  Company 

5.  North  American  Aviation,  Inc. 

6.  Northrup  Aircraft,  Inc. 
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7.  Rand  Corporation 

8.  U.  S.  Army  Research  and  Development 
Sub -office  (rocket),  Fort  Bliss,  Texas 

Little  work,  however,  appears  to  have  been  done  on  the  derivation  of 
theoretically  optimum  control  equations  or  operating  parameters.    For 
this  reason  the  writers  determined  to  examine  a  long-range  automatic 
navigation  system  now  under  development  at  the  Instrumentation  Labora- 
tory, M.I.T.,  as  a  means  for  exploring  theoretically  some  of  the  general 
problems  of  automatic  navigation.    The  Instrumentation  Laboratory  proj- 
ect, which  is  being  executed  under  USAF  Contract  W33-038ac-13969,  is 
designated  "An  Automatic  Navigation  System  -  Project  Febe,"  and  will  be 
the  subject  of  a  report  soon  to  be  published.    This  project  has  for  its  pri- 
mary purposes  to  determine: 

1.  the  feasibility  of  long-range  automatic  navigational 
guidance  of  bomber  airplanes 

2.  useful  design  parameters  for  a  serviceable  system 
for  military  use. 

Considerations  other  than  military  and  tactical  indicated  that  the  sys- 
tem should  employ  solar  tracking,  a  magnetic  azimuth  system,  and  a  con- 
stant ground  speed.    The  system  is  allowed  to  make  no  contact  with  the 
earth.  Altitude  is  determined  by  the  use  of  a  barometric  altimeter;  the 
only  inputs  to  the  system  are  celestial  observation,  the  magnetic  field  of 
the  earth,  and  the  various  accelerations  experienced  by  the  airplane. 

This  system,  which  has  been  installed  in  an  Air  Force  B-29  bomber, 
causes  the  airplane  to  fly  a  definitely  programmed  great  circle  course  at 
constant  ground  speed.   Vertical  accelerations  are  assumed  to  be  of  such 
small  importance  that  they  can  be  neglected.    Figure  1-2  shows  a  functional 
diagram  of  the  Febe  automatic  navigational  system.   During  the  period  from 
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June  to  August  of  1948,  the  writers  of  this  thesis  were  concerned  with  the 
Febe  project  in  an  under -instruction  status.    Through  the  courtesy  of  Dr. 
C.  S.  Draper,  they  were  also  permitted  to  attend  a  Seminar  on  Automatic 
(Celestial  and  Inertial)  Long-Range  Guidance  Systems  conducted  at  the 
Massachusetts  Institute  of  Technology  by  the  Scientific  Advisory  Board 
from  February  1  to  February  3,  1949. 

This  background  naturally  led  the  writers  to  concern  themselves  with 
the  form  which  a  practicable  service  guidance  system  might  take.    The 
Febe  system  presents  several  technical  and  military  problems,  which 
arise  principally  in  connection  with  the  following: 

1.  Constant  ground  speed 

2.  The  necessity  of  solar  tracking 

3.  Magnetic  azimuth  input 

4.  Weight  and  size  of  components. 

Since  a  successful  automatic  navigational  system  is  ultimately  des- 
tined to  fly  in  a  long-range  high-speed  missile,  it  must  be  assumed  that 
the  weapon  will  be  operating  at  close  to  the  maximum  range  permitted  by 
its  size  and  fuel  capacity.    The  navigational  system,  therefore,  should  not 
seriously  lower  the  fuel  economy.    This  suggests  a  system  which  will  fly 
at  constant,  or  nearly  constant,  airspeed.    Furthermore,  if  the  power  plant 
consists  of  an  athodyd,  the  airspeed  must  be  maintained  very  nearly  con- 
stant by  the  fundamental  limitations  of  this  type  of  propulsion. 

The  limitation  of  the  Febe  system  to  a  constant  ground  speed  was 
largely  dictated  by  the  requirements  of  the  azimuth  system.    The  magnetic 
input  could  be  eliminated  through  the  use  of  two  star-tracking  telescopes. 
The  problems  of  celestial  tracking  become  increasingly  severe,  however, 
as  missile  speed  is  increased.   At  supersonic  and  near-sonic  speeds,  it  is 
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no  longer  feasible  to  have  an  astrodome,  in  which  to  house  the  tracking 
unit,  projecting  from  the  fuselage.    Thus,  tracking  through  a  flat  window, 
with  the  attendant  difficulties  which  this  entails,  becomes  a  necessity. 
Also,  at  high  speeds,  thick  boundary  layers  and  intense  heating  will  exist 
along  the  surface  of  the  missile,  greatly  complicating  optical  problems. 
Finally,  the  inclusion  of  a  celestial  tracking  system  materially  increases 
the  space  and  weight  required  by  the  navigational  system  —  and  weight  and 
space  factors  become  increasingly  important  with  increasing  range  and 
bomb  load  of  an  aircraft.    It  appears,  however,  that  approximately  a  thirty  - 
fold  decrease  in  the  uncertainty  levels  of  existing  gyros  would  permit  the 
maintenance  of  an  inertial  coordinate  reference  system  within  the  missile 
through  the  use  of  such  gyros  alone,  thus  eliminating  the  use  of  celestial 
tracking. 

Great  progress  is  being  made  in  the  improvement  of  gyros.    The  work, 
for  example,  that  is  being  done  by  the  group  at  the  Massachusetts  Institute 
of  Technology  under  Dr.  C.  S.  Draper  shows  promise  of  obtaining  the  nec- 
essary accuracy  within  the  next  few  years. 

The  model  of  the  Febe  automatic  navigational  system  presently  in- 
stalled in  a  B-29  is  very  bulky,  and  weights  1917  lbs.    (See  Fig  1-3). 
It  practically  fills  the  after  pressurized  compartment  of  the  airplane.    Of 
course,  it  must  be  remembered  that  no  effort  was  made  to  decrease  the 
size  of  this  equipment,  or  to  keep  its  weight  at  a  minimum.    Nevertheless, 
a  missile  held  to  a  constant  airspeed,  rather  than  to  a  programmed  ground 
spped,  could,  by  using  gimbal  solvers,  greatly  simplify  the  elaborate  trigo- 
nometric and  speed  computers  which  Febe  requires.    Much  of  the  work  of 
the  computation  systems  might,  of  course,  be  performed  on  the  ground 
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prior  to  the  start  of  the  flight,  and  be  fed  to  the  system  through  tapes,  but 
this  would  introduce  further  complications. 

From  all  of  these  considerations,  a  self-contained  inertial -gravita- 
tional system  flying  a  great  circle  course  at  nearly  constant  airspeed 
appears  to  provide  one  of  the  most  satisfactory  solutions,  if  such  a  system 
can  be  implemented. 

After  this  examination  of  the  Febe  system,  there  remain,  then,  for 
use  in  the  purely  inertial  system,  a  gyro-stabilized  inertial  platform,  a 
clock  mechanism  to  remove  the  earth's  diurnal  rotation,  and  two  mutually 
perpendicular  single -degree -of -freedom  pendulous  accelerometer  units. 
In  addition,  there  are  gimbals  and  servomechanisms  isolating  the  inertial 
platform  from  the  motions  of  the  missile  containing  it,  and  orienting  the 
platform  so  that  its  axis  is  parallel  to  the  polar  axis  of  the  earth,  so  that 
a  reference  vertical  parallel  to  that  of  the  point  of  departure  (or  any  other 
convenient  reference  direction)  is  maintained  independent  of  the  motion  of 
the  missile.    (See  Figs    1-4,  1-5  and  1-6.) 

The  three  integrating  gyros  which  are  used  to  maintain  a  gimbal- 
mounted  inertial  platform  "fixed"  in  inertial  space,  parallel  to  the  earth's 
equatorial  plane,  are  mounted  so  that  they  detect  motion  of  the  platform 
with  respect  to  inertial  space  about  the  polar  axis  of  the  earth,  and  about 
two  mutually  perpendicular  axes  parallel  to  the  plane  of  the  equator  (Figs. 
1-7,  1-8).   Servomechanisms,  controlled  by  these  gyros,  rotate  the  gimbals 
so  that  the  inertial  platform  remains  fixed  in  inertial  space.    The  earth 
platform,  mounted  on  the  inertial  platform  with  the  identical  polar  axis,  is 
rotated  about  that  axis  once  in  twenty -four  hours,  to  transform  the  inertial 
reference  system  of  the  inertial  platform  into  an  earth  reference  system. 

8 


\ 


RANGE   INDICATION 
PENDULUM    UNIT 


CONTROLLED   MEMBER 
GIMBAL    SYSTEM    BASE 


TRACK    CONTROL 
PENDULUM     UNIT 


CONTROLLED   MEMBER 


CONTROLLED   MEMBER 
OUTER    AXIS   GIMBAL 


INNER 
ISOLATION 
G 


GIMBAL 


/ 

1 

0 

1        1 

3                 6 

INCHES 

1 
9 

1 
12 

FIG.   1-4. 
PHOTOGRAPH   OF    INERTIAL    REFERENCE    AUTOMATIC 
NAVIGATION    SYSTEM   MODEL. 


CO 


Q»R)(tcpl 


PAlfclep) 


^Z(GCHdep) 


'Elsetl 


QlIR) 


icmlitc) 


TRACK 
CONTROL 
SYSTEM 

PENDULOU' 
UNIT 


(S9)(DC)(ind)ltcU 
X 


TRACK    CONTROL 
MECHANIZATION 


'lcm| 


(5<))Az(cm)lind) 
5<fAz(cm)(true) 


COMPARATOR 


ISqJClA 


Z(cm) 


ROLL    AND  YAW 
DRIVE   SYSTEM 


COMPARATOR 


CD 


CONTROLLED 
MEMBER 


(cm) 


GIMBAL 
SYSTEM 
(FOR 
CONTROLLE 
MEMBER 


T 


lEPlles) 


EARTH 
PLATFORM 


JlPlllS) 


WE 


COMBIMOR 


(Sq)iCD) 


dndl 


t 


(Sg)(PAt 


_  tinqi 


EARTH  POLAR 
AXIS  CONTROL 

DIRECTION 
ANOLE  RESOLVER 


cosA 


[BWt_-JCD|jiILd|L. 


CORIOLIS 

ACCELERATION 

CORRECTION 

SYSTEM 


(Sq)A 


rlmdl 


(cm|(ri| 


(npl 


RANGE 

INDICATION 

SYSTEM 

PENDULOUS 

UNIT 


(Sql(DC) 


(indlllonq^    LONGITUDINAL 
MECHANIZATION 


Ar 


Ildestl-ldepl 


TIME 


Ar 


(indl 


UPlUsi 


INERTlAL 
PLATFORM 


(DIP, 


(IP)  (is) 


=  P      -P    • 

1  (JPKIS)     M(LS) 


GIMBAL 

ISOLATION 

SYSTEM 

(FOR   INERTlAL 

PLATFORM) 

I 


M(is) 


GYRO 
SYSTEM 


(OP 


(1  Pill  SI 


MISSILE 

(WITH 
CONTROL 
SURFACES) 


_^6rt)(trueHtc) 


(Sq)Ar(md)  _# 


RANGE 
ANGLE 
DRIVE 
SYSTEM 


COMPARATOR 


^Mjyert)  jdestj - (V e r t Uind][ _ 

MEASURE    OF  DISTANCE  TO  DESTINATION 


SEE  GLOSSARY  (APPENDIX    i_)  FOR   DEFINITION  OF  TERMS 
AND  SYMBOLS  EXCEPT   FOR  THOSE  LISTED  BELOW. 

Q(IR)Ucp)     "INERTlAL  REACTION  OF  GRAVITY  IN  TRACK 
CONTROL  PLANE- 
VECTOR  POSITION  IN  ONE  PLANE. 

(l>\J(dep)  "POLAR  ANGLE  AT  POINT  OF  OEPARTURE. 

Az(GCXjJep)"AZ"v1u'rH  0F  ANGLE  OF  THE  GREAT  CIRCLE 
TRACK  AT  THE  POINT  OF  DEPARTURE 


>E(set) 


(Sq) 

(is) 
(es) 

W1E 

(D) 
(0 


RADIUS  OF  EARTH  SETTING. 

-DENOTES  A  SIGNAL 

-DENOTES  INERTlAL   SPACE 
-  DENOTES  EARTH   SPACE 

-ANGULAR   RATE  OF  ROTATION  OF  EARTH 
WITH  RESPECT    TO  INERTlAL   SPACE. 

-OEVIATION 
-CORRECTION 


INPUTS 


^-    SIGNAL   OUTPUTS 

»       PHYSICAL    CONNECTION 
^~    POWER    LEVEL    OUTPUTS 


FIG.I-5.  SIMPLIFIED    FUNCTIONAL  DIAGRAM   OF  AUTOMATIC  INERTlAL  GUIDANCE   SYSTEM. 


10 


CONFIDENTIAL 


(PA) 


EARTH    POLAR  AXIS 


(dep) 


GREAT   CIRCLE   COURSE 
BETWEEN     DEPARTURE 
ANO    DESTINATION 
POINTS 

DESTINATION 
POINT 


CONTROLLED 
MEMBER 


DEPARTURE 
POINT 

(d«P>  /      AXIS   OF 

NCC    VERTICAL)  /OUTER  GlMBAL 

OP  CONTROLLED 
MEMBER 


FIGURE    1-6     GEOMETRIC    DIAGRAM   OF  AXES  OF 
AUTOMATIC    NAVIGATION    SYSTEM 


u 


SINGLE -DEGREE -OF - 
FREEDOM   GYRO  UNITS 


NERTIAL    PLATFORM 


OUTER 

ISOLAT  ION 

GIMBAL 


BASE 


FIG.    1-7. 

PHOTOGRAPHIC   VIEW   OF    INERTIAL    REFERENCE    AUTOMATIC 

NAVIGATION  SYSTEM  MODEL   WITH    POLAR    AXIS  TIPPED  AWAY 

FROM   THE    NORMAL    DIRECTION    TO    THE    BASE. 


ISOLATION    CONTROLLED   MEMBER 


/                                                                        -; 

r 

0 

1 

3 

1 
6 

INCHES 

1 
9 

1 
12 

N   GIMBAL 


BASE 


FIG.  1-8. 
PHOTOGRAPHIC  TOP   VIEW   OF   INERTIAL   REFERENCE    AUTOMATIC 
NAVIGATION  SYSTEM   MODEL   WITH   CONTROLLED  MEMBER  GIMBAL 

SYSTEM    REMOVED. 


CGNFJDEN  HAL 

Upon  the  earth  platform,  in  turn,  the  controlled  member  is  mounted 
in  gimbals  (Fig   1-9).    The  axis  of  the  outer  gimbal  of  this  member  is 
made  to  lie  in  the  plane  of  the  great  circle  connecting  the  points  of  de- 
parture and  destination  on  the  surface  of  the  earth,  considered  as  a 
sphere.*    Then,  about  the  inner  controlled-member  gimbal  axis,  the  con- 
trolled-member  is  moved  until  it  lies  parallel  to  the  surface  of  the  earth 
at  the  point  of  departure.    The  inner  controlled-member  axis  is  placed 
perpendicular  to  the  programmed  great  circle  path  so  that,  as  the  missile 
moves  along  the  programmed  path,  the  controlled  member,  by  rotation 
about  the  inner  axis  alone,  can  be  maintained  at  all  times  parallel  to  the 
surface  of  the  earth  below  the  missile.    The  angle  through  which  the  con- 
trolled member  has  rotated  will  then  be  a  direct  reading  of  the  distance 
travelled  by  the  missile  along  the  programmed  great  circle.   At  some  pre- 
determined angle  between  the  controlled  member  and  the  gimbal  connect- 
ing it  to  the  earth  platform,  the  missile  may  be  said  to  have  arrived  at  its 
destination. 

The  two  pendulous  units  are  mounted  upon  the  controlled  member,  one 
with  its  plane  of  motion  in  that  of  the  programmed  great  circle,  the  other 
with  its  plane  of  motion  perpendicular  to  that  great  circle  plane.    The  lon- 
gitudinal pendulum,  i.e.  the  unit  with  its  input  plane  parallel  to  the  plane 
of  the  programmed  track,  sensing  longitudinal  accelerations,  produces  an 
output  signal  which  is  then  modified  through  the  implementation  of  the  lon- 
gitudinal mechanisation  equation.   The  resultant  signal  is  used  as  the  input 


*  This  orientation  is  accomplished  by  rotating  the  controlled  member  on 
the  earth  platform  to  the  correct  meridian  angle  of  the  great  circle  track 
at  the  point  of  departure,  LAz(o,c.)(dep)^  and  by  mtin%  ^e  outer  gimbal 
axis  to  the  correct  polar  angle  (PA)/.onv  of  the  great  circle  track  plane 
with  respect  to  the  polar  axis  of  the  *aep;  earth. 
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to  a  servodrive  system  that  rotates  the  controlled  member  so  that  It 
tracks  the  local  vertical.   Similarly,  the  track  control  pendulum,  sensing 
lateral  accelerations;  (i.e.  accelerations  perpendicular  to  the  great  circle 
track  plane),  produces  an  output  which  is  modified  through  the  implemen- 
tation of  the  lateral    mechanization  equation,  and  the  resultant  is  used  to 
position  the  missile  control  surfaces  so  that  the  missile  remains  close  to 
the  programmed  great  circle  track.    The  altitude  of  the  missile  above  the 
surface  of  the  earth  is  maintained  by  a  pressure  type  altimeter  that  sup- 
plies the  essential  input  to  the  altitude  control  system. 

A  Coriolis  acceleration  computer  is  used  to  compensate  for  the  ef- 
fects of  the  Coriolis  acceleration  upon  the  action  of  the  system.    Geodesic 
acceleration  is  in  general  small  enough  to  be  ignored.   It  will  be  the  pur- 
pose of  this  thesis  to  discover  whether  it  is  theoretically  possible  to  main- 
tain the  errors  between  the  indicated  and  actual  local  verticals  small 
enough  to  permit  sufficiently  accurate  indication  of  the  local  vertical  for 
the  long-range  bombing  and  missile  guidance  problems.   Various  mechan- 
ization equations  will  be  examined,  to  see  which  offer  the  greatest  chances 
of  success.    The  most  promising  of  these  will  be  furthur  examined  to  de- 
duce the  optimum  control  parameters. 
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CHAPTER  H 
KINEMATICS 

1 .    Purpose 

It  is  the  purpose  of  this  chapter  to  derive  the  fundamental  kinematic 
equations  describing  the  motions  of  the  missile  postulated  in  Chapter  I.    The 
principle  geometrical  relationships  are  illustrated  in  Fig  H-l,  and  are  shown 
symbolically  in  Fig  n-2. 

The  following  simplifying  assumptions  are  made  in  deriving  the  kine- 
matic equations: 

a.  Lateral  and  longitudinal  motions  decoupled 

b.  Constant  altitude  flight  path 

c.  Spherical  and  homogeneous  earth 

d.  Controlled  member  located  at  center  of  gravity 
of  missile 

e.  Zero  angle  of  attack  and  side  slip  angle 

The  decoupling  of  the  lateral  and  longitudinal  motions  divides  the  prob- 
lem into  one  of  control  of  the  missile  in  the  direction  perpendicular  to  the 
programmed  great  circle  and  one  of  indication  of  the  distance  travelled  in 
the  programmed  great  circle.    The  kinematic  equations  are  therefore  de- 
rived separately  for  lateral  and  longitudinal  motions. 

Since  the  flight  path  is  to  be  considered  as  having  constant  altitude  above 
the  surface  of  the  earth,  the  only  vertical  component  of  acceleration  is  the  in- 
ertia reaction  acceleration  of  gravity.  With  the  earth  assumed  to  be  spherical 
and  homogeneous,  this  component  is  a  constant  over  the  surface  of  the  earth, 
with  the  numerical  value  of  115,920  feet  per  minute  per  minute  (equivalent  to 
32.2  ft/ sec  ).  The  assumption  of  a  spherical  earth  also  eliminates  any  geo- 
desic acceleration  terms  from  the  kinematic  equations  (ref  par  1  App  A). 
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TRUE  VERTICAL 


V«r4)(dvn) 


l(md) 

INDICATED  VERTICAL 


C)(V«HJ 

INDICATED    VERTICAL   CORRECTION 
DYNAMIC   VERTICAL  CORS 

(DC)(,nd) 

INDICATED    DYNAMIC   CORRECTION 


9ft  e) 


INERTIA    REACTION 
ACCELERATION 

JF    GRAVITY 


°ENDULUM 


COMTK'JLLED  MEMBER 


I  (c  my  vert}  — ' 


DYNAMIC    VERTICAL 


APPARENT     VERTICAL 
(PARALLEL  TO  DIRECTION  INDICATED  BY  PENDULU^" 

-(DC)^rue) 

TRUE    DVNAMK,    CORRECTION 


VERTICAL    ACCELERATION 
Cf   CONTROLLED   MEMBER 


u 


Ql  (ombres) 

8ESULTANT  ACCELERATION  OF 
CONTROLLED  MEMBER  WITH  RE- 
SPECT   TO    INERTiAL    SPACE.. 


t(crrtK  rescind  1 

NDICATED    RESULTANT 
ACCELERATION     OF    CON- 
TROLLED  MEMBER 

~Ql(cm)(hor) 

HORIZONTAL     ACCELERATION 
CF    LCNTROLLED  M£M3Ee 

<*Uc~; 

ACCELERATION   OF  CONTROLLED    MEMbtR 
WITH   REdPECT    TO   INERTIAL    SPACE  ;  ■',  IM- 
POSED   OF    ACCELERATION   OP     CONTROLLED 
MEMBER    WITH    RESPECT    TO   EARTH  {  N  - 
^LUDINO      GEODESIC,   e    HC  R  _:nta\_     CUCVA- 
~JRE,  ACCELERATION)  AND  ACE'-ERAT^N    OF 
COR. IO  LI  5 


•Q(Cor)*6(GEO) 

ACCELERATION   OF    CO<iOLi_> 
PLUS   GEODESIC    ACCELERATION 


TR'JE    vrv-TlCAL -DIRECTION)  OF  INERTIA  REACTION  ACCELERATION  OF  GRAVITY  AT  CONTROLLi  .    MEME"^. 

DYNAMIC    VERTICAL- DIRECTION  OF  VECTOR    RESULTANT  OF  INERTIA  REACTION    ACCELERATION  OF  GRAVITY, 

ACCELERATION    OF  COR.IOLIS    AND    GEODESIC    ACCELERATION    AT    CONTR;^_EC  'MEMBER.. 

APPARENT    vERTICAL-OlRECTION    OF   RESULTANT    ACCELERATION    OF    CONTROLLED    MEMBER     WITH   RESPECT 

TO   INERTiAL     SPACEC 

INDICATED  VERTICAL-DIRECTION    FIXED    TO  CONTROLLED    MEMfoEk.  -  PAR  ALLEL     TO    DIRECTION   OF    PENDULOUS 

ELEMENT      WHEN     ITS    OUTPUT     SIGNAL     I?    ZExO    WITH    CONTROLLED    MEMBER 
STATIONARY. 


FiG.n-a. 

GEOMETRICAL  RELATIONSHIPS   AMONG  DIRECTIONS   ASSOCIATED  WITH  THE  DIFFERENCES 
BETWEEN    THE   INERTIA    REACTION   ACCELERATION    OF    GRAVITY    AND    THE  DIRECTION 

OF  ACCELERATIONS    ON  A  MOVING   BASE. 
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A  Coriolis  acceleration  computer  is  included  in  the  fundamental  guidance 
system  (Fig  I~5),  and  is  assumed  to  work  ideally  (ref  par  2  App    A      ).    The 
kinematic  equations,  therefore,  do  not  contain  the  Coriolis  acceleration. 

With  the  controlled  member  located  at  the  center  of  gravity  of  the  mis- 
sile, the  gimbal  isolation  system  can  be  considered  to  remove  entirely  all 
effects  of  the  roll  and  pitch  of  the  missile.    No  accelerations  reach  the  con- 
trolled member  as  a  result  of  angular  accelerations  of  the  missile  about  its 
body  axes. 

Considering  both  the  angle  of  attack  and  the  angle  of  side  slip  as  zero, 
the  velocity  vector  of  the  missile  always  lies  along  the  longitudinal  body 
axis  of  the  missile;  that  is,  Vj/  .  ^  _  Mi  lies  along  Xj-  (see  Fig  n-3).    This 
is  taken  as  the  control  direction  (CD).    The  kinematic  equations  for  lateral 
motion  are  derived  first. 

2.    Kinematic  Equations  of  Motion  Perpendicular  to  Great  Circle  Track* 

The  geometrical  elements  of  the  simplified  lateral  guidance  problem  are 

illustrated  in  Fig  n-4.    Figure  II-5  defines  the  directions  and  angles  that  are 

used  in  the  derivation  of  the  kinematic  equations. 

As  can  be  seen  from  Fig  II-5,  the  following  relations  are  true: 

(C^A(vCrtxtc)  =  (DC)ctrueUtc)- (DC)((_ncjHtc)  (n-i) 


ton(DC)(traeKtc)  = g^ (n-3) 

♦This  development  follows  that  of  John  Hutztnlaub,  Notebook,  July,  1947. 
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CORRECT 
POSITION 


Z(cm) 


(C)[(Vert)(opp)(tc)3  cindHDQcindKtc) 
(opp)rtc)l(trucf^\true)(tc) 


FIGURE  IE- 5.     TRACK  CONTROL  DIAGRAM 

2^ 


^lE-(cm)Htc)  *  VLE_Cair)1(tc)  -  VL(aiP).M]sinA2.(crn)  cos(C)A(Vert)(tc) 


(n-4) 


From  (n-2)  and  (n-3) 

RE(C)A(Vert)(tc)    r  -9lRtQn(DC)(true)Ctc) 


(H-5) 


But,  from  (II- 1) 

tan(DC)UrueKU) 


=  tan[(aA(VerlHtc)    +  (DC)(indKtc)] 

tan(C)A(VcrtKto     +  tan  (DC)(inch(tc\ 
1-  tan  (C)  A(Vert)(to     to-ri^DOandiao 


(n-6) 


Therefore, 

P 

<C)A 


(VertXtC) 


IIR 


ton(C)ACVert)(tc) 


l-tan(C)A(VertUtc)     tan  (DC)(Cn(j)(tc) 


tan(DC)(indKt 


etc) 


l-tan(C)ACVertHtc)     tan(DC)(tnd)(tc) 


(H-7) 


This  is  a  theoretically  correct  expression,  if  no  vertical  accelerations  are 
present,  or  if  vertical  accelerations  are  assumed  to  cause  a  variation  in 
gjr,.    However,  if  the  size  of  (C)A/ ver+\/* c\  does  not  become  greater  than  10 
milliradians,  an  error  of  no  more  than  one  part  in  30,000  is  introduced  when 
(C)A/       .v  is  substituted  for  tan(C)A/ vertw+ cy    This  should  be  true  for  all 
cases  of  practical  interest.    If  this  is  done,  eq  (II-7)  becomes: 


34 


c°WH>EJVr/AL 

ir\A  .        /9iR   \  (C>AlVgrtXtc) 


nd)(tc> 


/9ir\     tQn(PC)anduta 

~^RE  '     l-(C)A(VertK1fc)    tan  (DC) 


(uxi)(tc) 


(n-8) 


If  -s~-   is  defined  as  W^,,  then  eq  (n-8)  becomes: 
E 

/rvA  +  W*    (C)A(VertUt0 

IWMW        -1   l-CC>A(Vert)rtc,  ton(DC)(indKtc( 


w*    tQn(DC)(lndutc) 

"E    l-^A(Vert>(tc)    ton(DC)(tndutc) 


(n-9) 


If  (DChindHtc)  *s  sma^er  tnan  f^teen  degrees,*  (DC)/.*!)/*-)  can  ** 
substituted  for  ^(DC)/.^/^  for  first  order  theoretical  solutions.   Also, 

(C)A/ vert\/f c\  *an(^-')(indHtc) '  as  ***  Proc*uct  °*  two  small  angles,  becomes 
very  much  less  than  unity,  and  can  be  neglected.   Then: 


2 


(C)A(Vert)(tc)   ♦  WnE  (C)A(VertKtc)  =  ~WnE(DC)(ind)(tc)  m-10) 


Using  the  same  small  angle  assumptions,  and  the  additional  assumption 
that  A„/ cm\  remains  less  than  fifteen  degrees,  eq  (n-4)  can  be  written  in 
the  approximate  form: 

V[E-<cm)](tc)   *      [E-(air>](tc>    ""  V[(air)-M]  Ai(cm) 

(n-ii) 


*  Ref  par  3  App  A 
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Equations  (n-2),  (n-10)  and  (n-11)  are  the  basic  kinematic  equations  that 
are  used  in  all  the  derivations  of  Chapter  in. 

3.    Longitudinal  Motion  Kinematic  Equations* 

The  longitudinal  problem  is  one  of  range  Indication.    The  geometrical 
elements  of  the  simplified  range  indication  problem  are  illustrated  in  Fig 
II-6.    Figure  n-7  defines  the  directions  and  angles  that  are  used  in  the  der- 
ivation of  the  kinematic  equations. 

As  can  be  seen  from  Fig  n-7,  the  following  relations  are  true: 

^DC^  (trueKlonq)  =     (C)A(Vert)(!ong)  +  lDC)Cind)(long)  (II"12) 

a 

'r(true)   ~ 


'LE-(cnrO](lon<j) 

A, 


(n-13) 


tan(DC)(trueHlo     f 1_ 

TIR 


Taking  the  tangent  of  eq  (n-12) 


(n-14) 


+~    <r\r\  tan(C)ACvert)(lonq)"»'ton(DC)CLadUlona) 

Tan  I UU  \true)i  lonq) — 7~T "XTT 3 ~ r=r^\ 

l-tan(C)A(Vert)(|oncp  tanCDC)anciK|o«S)  (n-15) 

Equation  (11-15)  is  a  theoretically  correct  expression  but,  as  in  the 
lateral  case,  all  practical  conditions  find  (C)A/ ver*\n onff\  smaller  than  ten 
milliradians.    Then,  eq  (n-15)  simplifies  essentially  to: 


♦This  development  follows  that  of  John  Hutzenlaub,  Notebook,  July,  1947. 
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(G|(Vert)(ind>(lonql| 

=  (C)  A(vert)(|0nqf(C)Ar 


GREAT 
CIRCLE  TRACK 


L"QE(crr7)fre5)(lonq)J 

(C)[CVert)(Qpp)(lori<3)J  (lnd) 


CVert)^d)(loncj) 

^   (cry,) 


I(cmKres)Uong) 


[E-(cml]  Oong) 


^[(Vert)(vnd)(\onS)]  (dep) 

=  A 


Wert)(dep) 

2 


^er^'f\rue)(\onq) 


[(Vert)(irueXlonc^]   Cdep) 

-  A 


fis.  n-i 

LONGITUDIMAL    INDICATION    DIAGRAM 
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C  0  N  i"  i  f) ! 
rwrv^  (C)A(VertH|onq^  tan(DC)andH|onq> 

WUfc'(trtttHlwq)        IWfTA  i  tan  (DC^ T~ i 


(H-16) 


From  eqs  (n-13)  and  (n-14), 

■•  •• 

+^%lT\n\  -       ^r(true)^E.  Ar  (true.) 

^  tir  wnE  (n-17) 


Substituting  eq  (n-17)  Into  eq  (H-16) 

Af(true)      _       CC)A(Vertu|onq)  +  ton(PC)(in<j)(lonq) 

W*nE  ^*^C)A(Vcrt)(|oa<1)   tan  (DC)(Lnju(on<5)  (U-18) 


From  this  there  results,  as  the  general  longitudinal  kinematic  equation, 
V  _  yy*   (C)A(VertK|onq) 

*r(true>  nE  f-OA^tHlon,)  ta"  ^C)andu,ons) 

^yy*    ton  (PC)and)(lon9> 

l-(C)A(Vcrtulon£^  tan(DC)and)(|onq)  (n-19) 


'nE. 


When  the  same  small  angle  assumptions  are  made  for  the  longitudinal 
kinematic  relationship  that  were  made  in  writing  eq  (n-10)  in  the  lateral  case, 
eq  (n-19)  reduces  to 


\(true)-  WnE  (C)A(VertUjon<^+  Wn£(DC)(tnd>(lon<j> 


(n-20) 
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This  equation  and  the  relationships 


(n-21) 


and 

Ar(app)   *  Ar(Lnd)  +  CDC)  CindKloncp  (11-22) 

form  the  basic  kinematic  equations  that  are  used  in  all  the  derivations  of 
Chapter  IV. 
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CHAPTER  III 
THE  TRACK  CONTROL  PROBLEM 

1.   Introduction 

The  simplified  kinematic  equations  developed  in  Chapter  II 
for  the  track  control  problem  are  restated  here  for  convenience. 


Ar«.  _<>  w,m  (+r\  — "  I^e  lw  A 


Cfc-(crnfl  (tc)        *»e  *W^(Vert)  ftc) 


(m-i) 


(QA(S/€rtXtc)  +Wut    (CJA^eHrXtc)-""*     Wa.    (DC)6„d)(tc)  (IH-2) 


*CE-(cm)](tc)  '    "  Vts-fdirflftc)       "Vcrair)-M\3     A 


item) 


(in-3) 


Using  these  simplified  kinematic  equations,  various  mech- 
anisation equations  are  investigated  in  an  effort  to  determine  a 
physically  realisable  and  reasonably  accurate  mechanisation  of 
the  track  control  problem.  As  previously  stated,  for  the  purposes 
of  this  study,  the  track  control  problem  is  considered  as  de- 
coupled from  the  longitudinal  problemfl)^  and  vertical  accelera- 
tions are  neglected!2).  It  is  further  assumed  that  the  missile 
operates  with  perfect  control  surface  servomechaniSms  and  that 


1  -  ref  par  2  App  A 

2  -  ref  par  4  App  A 
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it  possesses  perfect  aerodynamic  response^  '.    Finally,  the  pendulous 
accelerometers  are  considered  to  indicate  instantaneously  the  direction 
of  the  resultant  acceleration v  ' . 

The  purpose  of  the  track  control  mechanization  is  to  reduce  toward 

zero  any  angle  (C)  [(Vert)(mcn}(tc)  ^a*  nas  l)een  introduced  by  lateral 
motion  of  the  missile  away  from  the  prescribed  great  circle  track,  and 
to  maintain  this  angle  sufficiently  small  to  insure  arrival  within  one 
mile  of  the  destination.    The  control  is  achieved  through  operations  upon 
signals  received  from  the  measurement  of  the  available  angle 
(C)  [(Vert)app]/..x,  the  pendulum  angle,  in  order  to  correct  the  head- 
tag  angle,  A[X(cm)  -  (CD)]" 

Throughout  the  derivations  which  follow,  it  is  convenient  to  refer 
to  Fig  II-5,  which  indicates  the  relationships  among  the  important  angles. 

2.    Proportional  Control 

The  simplest  possible  mechanization  is  proportional  control,  in 
which  the  heading  of  the  missile  is  changed  in  linear  relation  to  the 
pendulum  angle  as  measured  between  the  pendulum  and  the  controlled 
member.    This  case  is  therefore  considered  first.    The  mechanization 
equation  for  this  concept  can  be  expressed  as: 

Azccnv*  =  ~  SLC(C0)][AA]  (DC)(Lndxtc)  (m-4) 


(3)  ref  par  5  App  A 

(4)  ref  par  6  App  A 
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Using  this  mechanization  equation,  and  the  simplified  kinematic 
equations,  numbered  one  through  three,  a  performance  equation 
for  this  system  is  now  derived. 

From  eq  (III-3), 


A-   ME-foir;3ftc)     VcE-fcm)](tc) 
ifcmF  w  (m"5) 

MTf'iilO-M) 


Integrating  eq   (III-l),   assuming  the  constant  of  integration  as 
zero^,   and  substituting  into  eq   (III-5), 


MJ 


From  eq   (III-2), 


HbCWr  fC^rKtc)  -HC)A(vJUtc)  (m-7) 


1  -  par  7  App  A 


S3 


LUNHD;:  ,  i  •- 


Substituting   (Ill-b)   and   (III-7)    into  the  mechanization  eq   (III-4), 


Vz.E-(&>r)mc)       *-£  V~'r*(sjert)(tc)    _ 


*  Tfair)- 


=  5 


CcfcDWCAA) 


[fair)-M3 


W 


(Vert)(tc) 


(iii-b) 


From  this,  the  performance  equation  can  be  written: 


E.W. 


E     »  'NE 


V        5 


^)Arvert)CtC)+  WNe  (C)A(vert)(tC)—  ~ 


LE-(&ir)]ftc) 


Kairi-M]  ^[c(cd)]Caa3 


TcCcdHCaA}        V£foir}-M] 


VT(Zir) 


■9) 


Although  this  equation  provides  a  stable  control  system,  a 
forced  error  results  from  the  lateral  component  of  the  velocity 
of  the  wind,  V r£  _  (air)]    »•  This,  obviously,  is  unsatisfactory 

for  the  control  of  guided  missiles,  because  no  large  steady-state 
errors  can  be  tolerated. 

3.   Integral  Control 

Next,  let  us  examine  a  system  that  changes  the  heading  of 
the  missile  in  accordance  with  the  integral  of  the  pendulum  angle. 
The  equation  for  this  mechanization  can  be  written 


Aifcm)"     5ccfcD):i[>^  /(DC)0„dK+c)  dt 


(111-10) 


n 


^Ohr  lOi  NiiAi 


L 


Differentiating  (III-6),  and  substituting  with  (III-7)  into 
(111-10)  differentiated, 

VC6-fair)]ftc)   »^E  (wAft/erQfo:)  _  c  l(CJA(ygrt)ftc)  .  ^.v  , 


(111-11) 


or 

(r\k         -I -^ fr)A  = — ^ \/ 

^'M(Vert)(tc)  '  a  ^W/VVert)(tc)     QIC         W        V[E-^r)]('tc)    (|||-12) 


.A]  Vc^ir).,, 


+  l 


This  system,  although  the  forced  error  is  now  due  to  the 
acceleration  of  the  wind,  has  no  damping  term.  S[(c)(CD)l  [AAJ 
can  be  made  large  to  reduce  the  size  of  the  forced  error,  but 
this  results  in  an  undamped  period  of  oscillation  of  approxi- 
mately 84  minutes  for  the  controlled  member.  Any  residual  forced 
error  (since  it  is  physically  impossible  to  make  S("/qwqD)1  rAjn 

infinite),  causes  this  oscillation  to  appear. 

4.  Second  Integral  Control 

In  order  to  reduce  the  forced  error,  consider  next  a  system 
that  changes  the  missile  heading  in  obedience  to  the  second  inte- 
gral of  the  penduluia  angle.  Here,  the  mechanization  equation  is: 


'Z(cm' 


^-•^co^  ff(DC\^Htc)  dt  dt 


011-13) 


d.Z 


CONI 

Differentiating   (III-6)   twice,    and   substituting  with   (III-7) 
into    (111-13),   differentiated  twice, 


Vfc -foiWtri        RE  (C)  f\(s,ert)kc)   _c  J(C/A^ert)(tc)    |     (r\*  \  CI    ,_14) 

(Yfcir)-M]  \.         WUE 


or, 


ff^A  -US  VcCd.r)-M3      />W  ,         <■       -  V[L<Air)-M3    /fU  L     \/ 


ie  '^6  ^e 


By  Routh's  stability  criteria^),  this  is  unconditionally- 
unstable. 

5.   Derivative  Control 

In  order  to  give  a  complete  presentation  of  possible  types 
of  mechanization,  consider  the  heading  angle  to  be  controlled  in 
accordance  with  the  derivative  of  the  pendulum  angle: 


A*,c*o  dt  =  -bLC(cvnck^    (DC)0Vw0ftc)  (iu-16) 


1  -  Routh,  "Advanced  Rigid  Dynamics". 
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Differentiating  (III-7),  and  substituting,  with  (III-6)  into 
(III-lo),  differentiated, 


coa 


(Vert^tc) 


^         ^(wA^)(tc)/*  (iu-17) 


From   (111-17) 
(C)A-      v   ^  +  W*   ll  H — l/hA  ~  Yli_^lM  ^ 

^^MvV.^c)T  wNE    I  -r  (UA^)uc)  - 

I     °fefco)]CAAa  Vrfair).MlJ  V[^)-Mj       v>rc.'c0)]lAA3 

By  Routh's  stability  criteria  for  a  cubic,  this  is  uncondi- 
tionally unstable.  Moreover,  there  is  a  forced  error  due  to  the 
velocity  of  the  wind. 

6.   Proportional  Plus  Integral  Control 

As  shown  in  paragraph  2,    proportional  control  produces  a 
solution  which,  although  it  has  a  forced  error  caused  by  wind 
velocity,  is  nevertheless  stable,  while  integral  control  (para- 
graph 3)  develops  a  solution  that  does  not  have  any  forced  error 
caused  by  wind,  but  which  also  has  no  damping,  and  is  therefore  on 
the  verge  of  instability.  Thus  it  might  be  expected  that  by  com- 
bining both  equations,  a  stable  system  without  forced  error  due  to 
velocity  of  wind  might  result.  The  resulting  mechanization  equation  is: 

™*<cm)~  -Wd)]Caa]  (^M-.'ndXtc)   ^CcCcd)]Iaa3  /  ( ^W(^d)^c)  Clt  (111-19) 
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Differentiating  eqs  (III-6)  and  (III-7)  and  substituting  with 
eqs  (III-7)  into  (111-19)  differentiated, 


•  •• 


=  5 


CO  A 


(\iert)(tc) 


[Cfco)3CAA3 


w: 


+(c)k 


Hi 


fvert)(tc) 


+  5 


;(oa 


cc(cd)j  taa:  J    \^j 


£^+(C)AfVt 


ME 


(V€rt)^c) 


(111-20) 


From  eq   (111-20), 


(C)A(^^hc)-\- 


vv  NE 


2 


CcfcDl](AAj        il    V^a.ir)-M] 


4- 


OccftoD 


Caa] 


w;E   j 


(C)Afvert)(*c)   +  WWE   (C)A(vert)^tc)  H~Wwe 


>5rcfcD)][AA] 


£(CD)]  [A  A]    V[&i7)-m] 


VrraiV)  -r 


V 


O  -  fiirfl        (111-21) 


By  Routh's   criteria   of  stability,    this  system  is  stable   if: 
to  StfofcoflCA*]    *M<i   5lwfcoarAA]    ^  pMitive  numbers, 


(b) 


w1 


~T<c) (c  oi]  Ca  a]    v  Ralr VM3 


>0. 


Furthermore,  as  expected,  this  system  has  no  forced  error 
resulting  from  wind  velocity,  although  it  has  one  caused  by  the 
acceleration  of  the  wind.  This  suggests  adding  second  integral 
to  the  system  just  discussed. 
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7.   Proportional  Plus  First  and  Second  Integral  Control 


Atfc^=-5fcfeDWx3CDC);ni)(te)—  SttfcBOtAfc/ftWtf«D(wd*     Sfcrcoi]DkXj //  (DC)rmd)wdt  dt 


Substituting  eqs  (III-6)  and  (III-7)  differentiated  twice, 
eq  (III-7)  differentiated,  and  eq  (III-7)  into  eq  (111-22) 
differentiated  twice, 


VL.«m-K.O)A 


fort)  fee)  


V 


~ ^[cfeomX  M   \         TT5-  ! "  lWMfo*)fcc) 


C(»'i>")-M] 


VV  WE 


-f-S 


®/W,  ,  fclv       1    5  fc)A 


fortKtc) 


i  w:, 


w 


NE 


roA 


(VertX-tc) 


(III -22) 


(111-23) 


From  eq   (111-23) 


9* 


4- 


S 


D^CA   A] 


S  V  s 

^Cc(co)3  {a  A3  vO»-)-nO         '-'ccfcoatA  A3  { 


(OK^r  w;«+ 


'CcIcoMCaaj 


S 

^tfcOUCA  A]   J 


(CXk  I       ^cfeP)JCAX3     vjl      /r^A  I      3ccfep)JCAA3    v^*     /£)a  =-- 


sc 


pfcpQCV  A3 


'Ccteo>3  Ca  aJ 


*'lcCcD>3  CAM      VCC&lr)-MT) 


VC£^,3       (HI  -24) 


According  to   Routh's  stability   criteria  for  a  quartic,   this 
system  is  stable  if: 

(a)   The  three  sensitivities  have  the  same   sign. 
(They  may  be  all  positive  or  all  negative), 


(b) 


W 


i. 

NE 


°Ccfc^CAA3        q  )  \  /  '      ^Cc(co)3t>A3  ^^ 


s     ....  V 

•^Lc-fco)  Laaj   I     Mj 


Laaj   I      tfa»*ii-)-M3 
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In  this  system  no  effect  of  the  wind  of  lower  order  than  the 
second  derivative  of  the  wind  velocity  enters  into  the  forced 
error.  It  seems  unlikely  that  a  forced  error  of  this  nature  could 
have  a  magnitude  large  enough  to  be  of  practical  importance.  How- 
ever, if  this  error  should  prove  important,  it  is  possible  to  go 
to  higher  order  performance  equations.  One  further  case  is  con- 
sidered. 


8.   Proportional  Plus  First,  Second  and  Third  Integral  Control 


Ancm)      ^ccfcouo; vn'DQ.'rvnfte)     S[cfcD)]Cx  a]  /( DC)r.nd)(tc)  dt    5Cc(co)K;^  / /TDQ,„dXte)  dtdt  5[cfcpaCA>;j // MDQ;h<j)fe)dt  dtdt 


Substituting  eqs  (III-6)and  (III-7)  differentiated  three 
times,  eq  (III-7)  differentiated  twice,  eq  (III-7)  differentiated, 
and  eq  (III-7)  into  eq  (111-25)  differentiated  three  times, 


(MI-25) 


V 


£(*ir)-M] 


w: 


fcplJCA  a'3 


1  (OAtwVtc) 


v^'^s/trdbdl     ^fcfcoa  ca  a] 


\@L^+(C)^:  H 


w1 

^  V*    WE 


VVe^Xtc)!        -TcfcDfl  TaaJ 


)(C)Arve^Vtc) 

w: 


<C)A 


/Vert)  (fa)/ 


(111-26) 


From  eq   (111-26) 


(C)  A^^fcsf 


S      V  . 


)-M3 


^-jos^  CC)ACve^-f  |w^  5c- 


♦JtcCcpOC' 


A    A} 


>j :  a  a'j 


"^fr'cDOCAADy 


fc)Arve^fe)4- 


Wi     OccCcpOCa  'A"J 


\(i 


'coi:  ca  a) 


5* 


—ass)  fc)AvtrtW+fw;  S 


CctcOlirA  Aj 


^tcfeolltt' 


sy 


•(C) A       +  {VnT     Ccfcp)3CAXj}/C)A        —  ~«   WNE  .-. 


V*  I  tj 


(IM-27: 


■"CcrooacA  AV 


S  \l  VI>-<W3 

C^coDC'X  Aj  VCTv'r)-Mj 
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By  Routh's  stability  criteria  for  a  quintic,  this  perform- 
ance equation  is  stable  if: 

(a)  All  coefficients  of  the  homogeneous 
equation  are  positive, 


(b) 


•TCtCPflrAA] 


w:- 


'CcCcdBC a  AD 
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[C(CD)DCAA1 


:c(co)3L'A'  a'J 
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f<3 
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(c(cp)3  L 'a  a3 


S        V  s 


IB. 


^CctoD)]  [AA3   *r(i!r) 
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^Cc(c 


O^tAA) 


M) 


'Ccko)]  CA  A3 


wil       j     ^cc(cp)3Ca  a 3    \ /\a/1        ^[c(cp)KaaJ      i       ^rcfco)DCAA] 


'NE 


'CC(CD)][A    A*J> 


N6   5 


[c(cd)DC'a  a  J 


'CCfcDDCA    A]> 


[cfcpUCA  A], 
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IR. 
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CcfcoHtA  A] 


5u      -..  V  S 

JMc(cdV)Ca  a  3   V[(»,r)-M3  ^Cc(cD)D[aa 


s  s 


'^(co)3C'a"a'J 
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C^coUCa  A'y 


'|R. 


S  V 

"^ Cc(c  013 La  'X]   v  K4 ir)- M] 


S      ..  V       s 

I  ^CctCPUlA  A3)     V/1    I     ~Ccdo)3CAA3 


5, 


[c(co)3L'a"a]' 


Ne  s 


5 


5 


Cc(cd)3Caa3 


Cc(co)3C'a  A3' 


E5  5 


IR 


Cc(CD)3tA'A'3 


L [Ql ' j 


M3 


'frfcp^ 


Tcfcp", 


These  restrictions  permit  the  sensitivities  to  be  either 
positive  or  negative,  provided  that  they  all  have  the  same  sign. 

It  is  of  interest  to  notice  that  in  each  stable  performance 
equation,  the  forced  error  has  the  same  type  of  coefficient,  with 
the  order  of  the  wind  derivative  increasing  from  velocity  for 
proportional  control  to  any  derivative  desired,  as  the  perform- 
ance equation  increases  in  complexity.  This  coefficient  is  of 
the  form: 


^CcfcPQEA  a  J      I       *^  rcfco)JCA  a'3  U    A    /  L    ^C=  (cd)2 C  a'a 


[c(cd«Ca  a  J 


'Ccfco)3[A  a] 


[cfepUCA  a ]/ 
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S      — V      s        \  /    s  s 


rv 
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S  V  S  /  V  v    5 


w 


E5  S 

Cc(cW3['a  a  J      ^fcteDXitTk'A'J    '«/ 


IR 


'frfcoQCA  A3 


5  V  S 


UCfCD^[*A"A3 


9 


IR 
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Cc(cp)3Ca  A3 


Cc(cp3]t'K  A]VCA 


1/  '     Q 

riAir)-M3  [cfcp)3UAi 


A], 


'CCfCPnCA  A] 


CCfCD)3tA  A3 
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The  forced  error  is  thus  seen  to  decrease  as  the  missile  airspeed 
increases,  and  as  the  sensitivity  S[c(CD)]  [M]  increases.  Theo- 
retically, if  this  sensitivity  increased  to  infinity,  there  would 
be  no  error,  either  transient  or  in  the  steady  state,  caused  by 
the  action  of  the  wind.  This  is  similar  to  the  conditions  found 

with  integral  control,  eq  (111-12)  whor*  a*  +h*  .    •«.•  • 

>  ^m  v-lxa  i.c  i ,  wnere,  as  the  sensitivity 

S[C(CD)]  [AA]  increases  to  infinity,  both  steady  state  and  trans- 
ient wind  errors  reduce  to  zero. 

^ — Track  Control  Correction  as  a  Function  of  Wi nH 

Since  the  mechanization  eq  (111-22)  for  proportional  plus 
first  and  second  integral  control  appears  to  provide  a  stable 
system  without  excessive  forced  error,  this  system  is  considered 
further.  The  performance  equation  for  this  system,  eq  (111-24), 
can  be  expressed  in  operational  form: 


'WA(v«rt)fte 


)  __ 


*  NE 


<S        w         r  (111-28) 


.  •-  7 ■ V  -tc(CP)JLA*]    »CCa.ir)-Ml     J 

vy ccccD)3 ca xj Vr&ir).M1     orwCDI7rvvJ'       (        s .....  I         <\  ne (r        C  W.. 


S       ....      NE  r       c 


,JLAAJ  W6' 


S> 


CO.  < 


If  (C)A/vertj  is  expressed  in  minutes  it  is  a  measure,  in 
nautical  miles,  of  the  transverse  linear  distance  of  the  missile 
from  the  desired  great  circle  track.  The  oresence  of  the  opera- 
tional symbol  in  the  numerator  of  the  equation  indicates  that 
there  is  no  error,  in  the  steady  state,  from  an  accelerating 
wind.  If  it  is  desired  to  find  the  missile  response,  A2(cm),  to 

an  accelerating  wind,  it  can  be  calculated  if  eq  (111-28)  is 

AZ(cm) 
multiplied  by  an  expression  for  -j~*J, ' The  result 

(c'A(vert)(tc) 

yields   the  heading  of   the  missile  as  a   function  of  wind    velocity. 
Manipulating   eqs    (III-l),    (III-2),    (III-3),    and    (111-22),    and 
expressing  the   result   in  operational  form, 
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Combining  eqs  (111-28)  and  (111-29)  gives, 
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The  negative  sign  shows  that  the  heading  Az(cm)  of  the  air- 
plane  is  correct  to  oppose  the  tendency  of  Vk  .  (air)]  *o   blow 
the  missile  off  the  great  circle  course.  This  is  true  even  if 
the  sensitivities  are  negative. 

10.  Track  Control  Closed  Loop  System 

The  relationship  of  eq  (111-30)  can  be  obtained  from  a  study 
of  the  over-all  lateral  system  as  a  servomechanism.  In  the  direct 
servo  closed-loop  system,  the  angle  (C)A/vert)  does  not  appear, 
but  it  can  easily  be  obtained,  as  indicated  in  the  accompanying 
Pig  III-l. 

In  this  figure: 

(FF)i  has  the  value 


SIR 


(PF)2  -    ziCM) and  is  made  up  of  the  sum 

*DC' (ind) (tc)   of  two  performance  functions, 


(pf).  -  .(^!llnd)(tc)     «* 

(DC>(true)(tc) 


(PF)D  m       Z[cm) which  is  the  mechanisation 

(DC) (ind)(tc)  equation  in  operational  form 


(PF)3,  in  the  feedback  path,  is  (-  Yt(air)  =  MJ)p> 

SIR 

providing  the  proper  function  for  the  com- 
binor,  which  is  physically  the  pendulous  unit. 
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The  pendulous  unit,  if  (DC)(true)  remains  smaller  than  about  fif- 
teen degrees,  solves  the  equation 

<3,R   '      ^    "^^"^  (,,,-3,) 

This  equation  is  derived  from  eqs  (III-l),  (III-2),  (III-3) 
differentiated,  and  the  relationship  that 

®oc±M= rc)A(vfert)fte)+  (oo(ini)M  cm  -32) 


This  last  equation  can  readily  be  obtained  from  an  inspection  of 
Fig  II-l. 

Solving  for  the  performance  function  (PF)a,  using  eqs  (III-l), 
(III-2),  (III-3),  and  (111-32), 

* 

The  mechanization  equation,  which  in  this  example  is  eq  (111-22) 
expressed  in  operational  form,  is 


s 

fPFl=   M^cm)    - 


Px4>-rfe£sagap4.  ^fcfcp)3I>*H  (111-34; 


tt^Ud)(ic)  P 
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The  product  of  eqs  (111-33)  and  (111-34)  gives 


(PF)  = 


A 


s    ...       s     l  r<3 


(DC) 


ftruueKtc) 


P1 


(in  -35) 


but,    from  the  theory  of   servoroechanisms, 


(1) 


A.,cw)    =  (PF),  fiPFj 

si*-**    i  +  rpFi  cpf)5 


(in-36) 


Solving  this  equation  gives  eq  (111-30),  as  before.  The  stability 

criteria  for  eq  (111-30)  are  identical  with  those  for  eq  (111-24). 

Since  the  sensitivities  can  have  either  sign,  a  necessary  condi- 

SIR 

V[E  -  (air)] 

if  the  sensitivities  are  chosen  as  negative.  If  this  is  expressed 

as 


tion  for   stability   is  that 


S[C(CD)][AA]        > 


5       ...  =  -  K 


'|R 


v: 


(wWeK>ia*<15<0,), 


£  -    air)] 


(Mi-37) 


1  -  Brown  and  Campbell,  "Principles  of  Servomechanisras". 
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then,  for  stability, 
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11.  Numerical  Solutions,  Using  Two  Quadratics 

Using  the  relationships  that  have  been  established,  and 
substituting  numerical  values  for  the  coefficients,  it  is  possi- 
ble to  solve  for  the  various  angles  as  a  function  of  the  accelera- 
tion  of  the  wind  Vrg   (air)]  •  From  a  sequence  of  such  solutions, 
using  a  range  of  values  for  the  sensitivities,  it  is  possible  to 
select  the  most  satisfactory  values,  and  to  estimate  the  effec- 
tiveness of  the  system.  This  is  now  done. 

The  homogeneous  portion  of  the  performance  equation, (111-24) , 
is  a  biquadratic.  It  can  then  be  written 


(px+  2(DR)W,  +W,1)  (p*+z(DRl  K  +Wtl)  -  0 


(111-39) 


which  factors  the  equation  into  two  quadratics.  Manipulation  of 

the  equation  indicates  that  if  the  damping  ratios  and  natural 

frequencies  are  varied,  the  higher  damping  term  associates  itself 

with  the  shorter  of  the  two  frequencies.  This  result  then  appears 

essentially  as  that  of  a  lightly  damped  long  period  quadratic, 

which  is  undesirable.  For  this  reason  (DR)^  and  (DR)2  ,  and  W^ 

and  W2  have  been  chosen  equal  for  all  solutions  of  the  equation, 
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and  equations  art  now  derived  to  facilitate  the  selection  of 
values  for  the  sensitivities  to  make  these  limitations  hold. 
Expanding  (111-39),  in  the  general  form,  gives 

p4+[2  (DR)W+Z  (DR\W]  p3+  [W,l+Wtl+  MO^D^y^p^ 

[2(D^yva+z(DR)xw;]p+w;w;.  =o 


(111-40) 
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(FR),    and    (DR)!   -    (DR)2   -    (DR)    then, 


comparing  the  coefficients  of  eq  (III-40)with  those  of  eq  (111-24) 
gives  the  following  relationships: 
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When  these  four  equations  are  solved  simultaneously  for  the 
sensitivities  and  for  the  damping  ratio,  the  results  are: 
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2(FR) 


(111-45) 


<3iR 

5wcoBC"J=  8WBEVtlM3CDR)x(FRr  ( " ' _46) 

* 

t£,"atiX!    2V0MPMflfl>R)  (,,,-47) 

a**,,*     8V^,.M](W  (l,|-48) 


It  is  noticed  that  the  choice  of  a  damping  ratio  immediately 
fixes  all  sensitivities  and  the  frequency  ratio  (or  that  the  choice 
of  a  frequency  ratio,  alternatively,  fixes  the  damping  ratio).  As 
the  frequency  ratio  increases;  i.e.,  as  the  quadratics  have  pro- 
gressively higher  natural  frequencies,  the  damping  ratio  simul- 
taneously increases  without  limit.  Conversely,  as  the  frequency 
of  the  quadratics  approaches  that  of  the  &4  minute  period  "earth 
pendulum",  WNE  ,  the  damping  reduces  to  zero.  The  sensitivities, 
at  the  same  time,  increase  without  limit. 
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12.  Numerical  Solution  Using  Four  First  Order  Terms 

A  second  method  of  writing  the  quartic  equation,  which  can 
be  examined  easily,  and  which  may  yield  results  of  practical 
interest,  is  that  in  which  the  equation  is  reduced  to  four  decay- 
ing exponentials.  This  equation  is  of  the  form 


(p+^,)(p+o)(p+a)(p+(F). 


—  0  (111-49) 


When  expanded,    this  becomes 


P^kcT)  Vr)/  (cT)3+fcT)Jp3+  [(cT)t(cr)+<cT\(cT\*  rcrycT)/  (ct\(ct)*  (ct)Jct\ 

"l"CCT)/CT)Jpa~h  [fCT\(Cr)t(CT\  +  (CT\(CT)X(CT)*  rCT)/CT)3  (CT^fCT^  (CT^CT)  J  P 

. ! 

(ct\(ct\(ct)3(ct)+  ~°  (111-50) 

The  general  relationships  among  the  coefficients  of  eq  (111-50) 
and  eq  (111-24)  give 
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9  S 


+ 


*S  V  s 

ycc  fco>]  CA  A3    v  Cfai  r)_  M J  "*  Cc  ( c  o  V}  z  a  a  3  _ 


(tD+CCtt+CcT),  -KC'TJL     ( m  -51 ) 


(MI-52) 
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(Hl-53) 


CCfCOWCA  Al 


"     (CT}(CT\(CT)JCTl 


(111-54) 


The  special  case  which  is  considered  here  is  that  for  which 
(CT)i  «  (CT>2  -  (CT)^  -  (CT)j^  -  (CT).  Under  these  special  circum- 
stances, the  simultaneous  solution  of  eqs  (111-51,  52,  53,  and  54) 
yields 

V^(G±V5z), 


(bf    z 


(in-55) 


from  which 


cct) 


1    =W„.  2.418  or^=W„  Q4IZ. 


'MC 


fcT) 


(111-56) 


ie 


CctcoXlCA  X3         W 

LfeW-M] 


w;. 


4w;    4 


CCT)     (CT): 


I 


I* 


CcCcoHCAA} 


V 


O^lD-M] 


(CT)X 


w;- 


52 


(CTl 


(m-57) 


(mi-58) 
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S  ,      -  = 


V 


C(a»»-)-Ml 


I 


(Of 


4W*— 


(at 


(IM-59) 


Evaluating  (111-57),  (III-5S),  and  (111-59),  using  eq  (111-56) 
with  the  smaller  (CT), 


IR. 


TcCcoUIaAi 


M^vM3WMt45 


(111-60) 


TcrcD))CAAD 


<3,R  1.205 
V 


( 1 1 1 -61 ) 


<3  W 


CcCco)DCaa3 


M,;^0728 


(III -62) 


This  gives  negative  values  to  the  sensitivities.  The  alterna- 
tive evaluation  would  make  the  sensitivities  positive,  but  would 
result  in  characteristic  times  fdr  longer  than  could  be  tolerated. 
The  system,  with  negative  sensitivities,  is  stable  by  Routhfs 
criteria. 
13 »  Numerical  Solution  Using  Quadratic  and  Two  First  Order  Terms 

The  final  method  of  examining  a  quart ic  equation  which  may  be 
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of  interest  is  that  in  which  the  quartic  is  broken  up  into  one 
quadratic  term,  and  two  decaying  exponentials, 

(p'+Z(DR)Wp  +  Wt)(p+(^:))(p+  g-j)  (in-63) 


Upon  expanding  eq  (III-63)  this  becomes, 

rgfoiflw  .     W.ji/'Ib  (111-64) 

lfcT),fcT^   fcT)/ fcir^Jr 

The  general  relationships  among  the  coefficients  of  equations 
(111-64)  and  (111-24)  are 

k H? +  |sfe£a£Aii]=  2cDR)w+-L.+J-)    , 

l^Kfc.ntAi]  VWr).M]  5ccfconr>iA]J  R     lCT7t  (111-65) 


5— 7w--^wn+(cT),+  ca)/ 


5 W"  =  <tf)(CT)  (m"68) 
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These  equations  can  be  solved  simultaneously  for  the  sen- 
sitivities, and  for  the  relation  which  connects  (CTJ^  ,  (CT)2 
W  and  (DR).  It  is  convenient  to  use  the  relationships:  (FR)  ■ 


W 


NE 


^  [(CT)PR] 


(CD^E 


2-ir 


and  [(CT)PRj2 


(CT)2WNE 


2-rr 


.  The 


resulting  equations  are: 


41T1  (FR)!|YCT)PR]  [(CT)PR]  +  I  =  4TT1  [(CV)  PR],  [(CT)PR\  +  4TT  (DR)(FR)  ([(CT)PR],  +[(CT)  PR],  }  +  I 


; 1 1  I -6?) 


ZTT1  Qx  fFKV  [fcT)PR],  [fcT)PR]t 


fc(CO>]  [A  A] 


Cc(co)]  CAa3 


H,lir).M]  W„£  (4TtI(DR)(FR)[(CT)PR],  [(CT)  PR\+ir(m  ([fCT)PR], +[CCT)PR]2  )-(DR)(fr)-TT  ([(CT)PRL],  +[fct)PR],')}    C 1 1 1  -7 

<3„  f(DR)(FR)-Hr([;(CT)PR]l  +  [fCDPE].  )} 

V.M]  [4TT2(DR)(FR)[(CT)PRj,  ^CDPRl  +  TffRy^CTjPiq+pCDPRlJ-COR^FIZ)-  IT  ([(CT)PR],  +  |fcT)PR|,j 


(M  1-71 ) 


<3«  w. 


rcfcDOtAX: 


V 


L(o»-nO 


(8trYDRXFR)LCCT)PK],  [fcflPRj,.  +  21T  fFE)1  ([(CT)PR],+RcT)PE]>2CDR)(FR)- 27T  ([fcT)PE],  +  [fCDPE],)];,,,^) 


CONFIDENT^ 


These  are  evaluated  for  the  single  case  of  (DR)  •  If   ,  with 
[(CTJPRJ!  -   [(CT)PR]2  -  0.1.  Then,  from  (111-69),  (FR)  -  0.307 
This  gives: 


S„X*,  y^  Wms    3764  (||1.73) 


9.J.205 


9,.  WHi  0.714 


^crcD^cAAj-  vi  (III  -75) 

^OairV-MiJ 


Notice  that  once  again  the  sensitivities  are  negative. 
14.  Plots  of  Results 

It  is  now  possible,  using  various  coefficients  for  the 
equations  of  paragraph  11,  or  the  indicated  coefficients  of  the 
equations  of  paragraphs  12  and  13 ,  to  obtain  numerical  answers 
by  solving  the  equations  with  various  wind  inputs*  This  was  done 
on  the  Rockefeller  Analyser  at  Massachusetts  Institute  of  Technology. 

It  is  not  necessary  to  solve  for  various  missile  speeds,  since 
the  speed  of  the  missile  can  be  combined  with  the  heading  angle  of 
the  missile  Az^caj#  This  is  possible  because  of  the  linearisation 
assumption  which  was  made  in  deriving  the  kinematic  equations,  and 
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holds  only  if  the  missile  airspeed  is  very  nuch  greater  than  the 
velocity  of  the  wind. 

Plots  were  made  of  the  angular  error  in  minutes  of  arc  of 
the  position  of  the  missile  perpendicular  to  the  great  circle 
track.  Considering  the  earth  to  be  spherical,  this  is  directly 
a  measure  of  distance  from  the  track  in  nautical  miles.  For 
convenience,  on  all  plots,  the  unit  of  angle  is  the  minute  of 
arc,  and  of  time,  the  minute  of  time.  The  foot  has  been  taken 
as  the  unit  of  distance* 

The  first  input  of  wind  which  was  used  in  solving  the  equa- 
tions was  a  rectangular  pulse  of  wind  acceleration,  or  a  constant 
acceleration  of  the  velocity  of  the  wind  from  aero  to  some  finite 
value.  The  duration  of  the  pulse  was  taken  as  one  minute.  Next, 
responses  were  taken  for  a  step  function  input  of  wind  accelera- 
tion. Finally,  responses  were  determined  for  sinusoidal  inputs  of 
wind  velocity  using  several  frequencies.  Plots  of  the  results 
follow,  using  various  sensitivities  in  eq  (111-22): 
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r^klL 


(DR) 

.417 

1 

uCctCD)HAA] 

Q 

°LC(CD)UAA1 

^[C(CO)]lAAl 
°[C(CD)][AA] 

^LC(.CD)HAA]    V[(air,-M] 

0.00000208 

0.4195 

0.0281 

.750 

0.00001158 

1.789 

0.0885 

1.050 

0.00003540 

4.890 

0.2172 

1.500 

0.0001255 

18.05 

0.6575 

1.875 

0.0002891 

35.75 

1.4210 

NUMERICAL  CONSTANTS  USED  IN  OBTAINING 
RESPONSE  OF  TRACK  CONTROL  SYSTEM  TO  PULSE  AND 
STEP  OF  WIND  ACCELERATION  AND  CIRCULATING  WIND 

TABLE  m-1 
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CHAPTER  IV 
THE  RANGE  INDICATION  PROBLEM 


1,   Introduction 

The  simplified  kinematic  equations  developed  in  Chapter  II 
for  the  longitudinal  problem  are  restated  here  for  convenience: 

vDwfinMlonaT  xJT.  ~~   (CjAjVertKlo^)  (IV-1) 


Arfiiwir  Ar(truO    I   (C)  A  r 


(IV-2) 


Arrapp)   ArM)+(DC;flnel)^ohJ) 


Using  these  simplified  kinematic  equations,  various  mech- 
anization equations  are  investigated  in  this  chapter  in  an  effort 
to  determine  a  physically  realisable  and  reasonably  accurate 
mechanisation  of  the  range  indication  problem.  The  assumptions 
are  identical  with  those  for  the  lateral  problem*  The  missile  Is 
assumed  to  move  with  instantaneous  response  to  wind  forcing 
functions. 

The  purpose  of  the  longitudinal  mechanisation  is  to  reduce 
toward  sero  any  angle  (C)  [(Vert) (ind) (long)]  which  has  been  in- 
troduced by  motion  of  the  missile  along  the  prescribed  great 
circle  track,  and  by  measuring  the  angle  between  the  reference 
vertical  and  the  indicated  vertical,  to  obtain  a  measure  of  the 
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location  of  the  missile.  The  orientation  of  the  controlled  member 
is  achieved  by  measuring  the  angle  (DC)  ( incj)  (long)  between  the 
position  of  the  pendulum  and  the  controlled  member,  and  by  oper- 
ating on  this  quantity  to  position  the  controlled  member  with 
respect  to  the  reference  vertical.  The  longitudinal  reference  is 
chosen  as  the  vertical  of  the  departure  point,  (Vert ) (dep) .  Since 
the  value  of  (DC)  ( incj )  (long)  is  dependent  upon  the  position  of 
the  controlled  member,  as  well  as  of  the  pendulum,  it  is  more 
convenient,  in  deriving  mechanizations,  to  use  the  angle  Ar(app) 
between  the  reference  vertical  and  the  indicated  vertical  of  the 
controlled  member,  than  to  use  (DC ) (ind ) (long)  directly.  Ar(app) 
is  physically  measurable,  since  it  is  the  sum  of  (DC)  ( j.nd)  ( long) 
and  Ar(j_n(j),  which  are  both  known.  It  has  the  advantage  of  being 
independent  of  the  position  of  the  controlled  member. 

Throughout  the  derivations  that  follow,  it  is  convenient  to 
refer  to  Figure  II-6,  which  indicates  the  relationships  among 
the  important  angles. 

2.   Trial  Mechanizations 

A.  The  first  mechanization  that  is  considered  is  that  in 
which  the  second  integral  of  Ar(app)  moves  the  controlled  member 
according  to  Ar(ind)  and  its  first  two  integrals.  That  any  simpler 
case  produces  an  unsatisfactory  solution  is  indicated  by  the  mech- 
anization equations  for  the  track  control  problem.  The  mechaniza- 
tion equation  is 


S[c«m>]CA«]  //  Ar(iBj)  dt  dt  +5[cfcml5 1«]  /  ARf;„j)  dt  +  ot 


ArrtnJ=  // Ar(apP,  dtdt 


[c(cw,YJ[Xa]    rVivwl  ~       //     Mr(app)    ULCJi»  (IV-4) 


..■i 


cor, 

■s 

Differentiating  twice, 

0[CfcmJ][AA]    r\r(\r\6)     '    ^[cCcw)]  [aaJ     ^r(md")      '    ^(cwiQCaaI     ^V(md)-  AK&pp)  (IV~5) 

Substituting   (IV-3)    into    (IV-4) 

Pccfcm)][AA]""  M    Arfind)    '     •^[c<'cmr)[AA3  AVindP"  ^[cfcm)}[AA'J     ArfiniJ)=  L/L(irs<j)|0^3 


(lv-6) 


which  is  the  mechanization  equation  in  terms  of  (DC) (ind) (lone) 
and  Ar(ind).  Further  substituting  (IV-1)  and  (IV-2)  into  (IV-6), 


/ 


l!  W*5 


I 


ME  ^[CCcmflCAA]  ; 


A  ^[crcm)J[AA3      A  Prc(cm)KAAT  '  I    A  _  Ma  I        ^Ce(c  rvO]  ZA  AT  (WAr  ^[cfcm)]  [AA]    (UAi 

ArCtrue)  -  rVftrue)         |  <.  c  Ar(trw)-^/nr  T  ~t" 

^LcrcvyilirAA]  l^[c(cmi:[AA)  ^[c(cwO][aa]J  ^[cfcm)]  [a  a]  0[c(c*vo;j 


This  control  equation  indicates  that  the  value  for  S  [c(CM)l  [  AAl 
should  be  unity,  and  it  will  be  so  chosen  in  further  considerations 
of  the  problem.  Even  with  the  position  error  eliminated  by  the 
proper  selection  of  S("(qwcm)1  [aa]  >  tnere  remains  an  error  caused 

by  the  rate  of  change  of  Ar(true).  Since  r(true)   ^3  a  measure 

"E 

of  the  distance  travelled  over  the  surface  of  the  earth,  and  because 
airspeed  is  constant, 


[AA] 

Ov-7) 


(IV-8) 


•  I 


.  < 


This  means  that  the  steady  state  error  is  a  function  of 
wind  velocity,  and  that  this  equation  is  therefore  unsatisfactory. 

B.   In  order  to  improve  this  response,  the  first  integral  of 
^r(app)  *s  added  to  the  second  integral  in  the  preceding  mech- 
anization equation.  The  resulting  equation  is 


A r(ind)dt<lt+ S[cfcmmAADy  Ar(imJ)dt  +-5Cc(on^A]  Ar„-nd)  -J   Arftpp)  dtdt  +  5Cc(cm,DCAA3  /  Arrap^dt 


(IV-;) 


from  which,  with  (IV-3), 


Further  substituting  (IV-1)  and  (IV-2)  into  (IV-10), 


(IV-10) 


J_    S^^.:   A-        _f|--L    _J U  =(C)Ar+|^^  (OAr  +^ fc)Ar 

l    ^  "rterue)     1     v./1    c         I  "r  true)  C  S  ,    ..  (.IV-llJ 


In  this  control  equation,  there  is  a  forced  error  resulting 

from  the  acceleration  of  the  wind  which  can  be  removed  if 

1  7 
n —  is  made  equal  to  W^E*  Under  this  condition,  a 

S[C(CM)][AAj 

stable   control   equation  results  with  a   forced   error  resulting 
from  no   effect    of    the  wind   of   lower   order   than   the   time   rate   of 
cnange   of  acceleration.    The  homogeneous   equation  contains  damping, 
but   has  a    period   of  about   34.4  minutes,    because   the   natural   fre- 
quency of  the   system  is  fixed  as  W2NE  when  the   value   of   S[C(CM}]  [AA] 


is    chosen.    Bee, use    it    is  desirable   to    increase   this   natural   fre- 
quency,   in   order   to   reduce   the   time   during   which   transient   errors 
are    of    importance,    other  mechanization   equations   are   examined    in 
this   chapter,    in  an   effort    to   acquire    control   over   period  as  well 
as   damping   in   the   control   equation. 

C.      The   third   mechanization  equation  to   be   considered   uses 
the   third    integral   of   Ar(app)    t0   move   the   controlled   member)    by 
operating  on   Ar(ind)    and    its   first    three   integrals.    This   mech- 
anization  can  be   expressed 


"rfINd)^^^%m)][w]    /  J^rtmd)  d"tdt  +  Orc(cm)J[KA3    /  Arfind)<n   H"   J[cfcrn)]r^A]  Ar(;nd) 


A,,app,Jtdtdt 


(IV-12) 


Differentiating  equation  (IV-12)  three  times,  and  substitut- 
ing equation  (IV-3)  into  ( IV-12], 


'[Cfc^flCAAJ 


^rCmcl)      :       ^> 


CC(C^)][AA] 


"i-fi'nd)     ""    ^Ccfcm)]CAA]    "rd'wd)    — ^^'(ind)(\o^ 


Further  substituting  equations    (IV-1)   and    (IV-2)    into   eq 
(IV-13), 


(iv-13) 


Vftrue) 


Cc(cm)]CAA] 


Tcfc^lCAAj 


I 


Wwe   0[cfcm>3[AA3 


I'lira.e) 


5, 


Cereal]  [A  A] 


^[C(C 


[c(cm)][AA] 


A. 


(irue) 


(C)Ar  + 


S* 


(CW)][AA] 


[C(cmVj[AA] 


(C)Ar  +    J[ot"*"i     (C)Ar  + 


(C)Ar 


'TcCcwOC  AA] 


5fc 


[cCcml][AA] 


This  control  equation  again  contains  a  forced  error  propor- 
tional to  the  velocity  of  the  wind,  and  when  the  proportionality 
factor  3r  ,-,][- AAj  i3  made  zero,  the  homogeneous  equation  becomes 


(IV-14) 


17 


NFIDENL'/.L 


unstable. 

D.  In  an  eifort  to  improve  this  situation,  the  second  in- 
tegral of  Ar(app)  is  added  to  the  third  integral  in  the  preced- 
ing mechanization.  Then, 


HI- 


A^cftdtdt-hS 


Rfcm)][AA] 


■/fa 


l<h 


dtdt+5 


Cc(cm)J  cXa^ 


I hr(\«t)  c"t+->cccc^]UAA]  Arf,-Wj)- -J  II  Ar(app)  dtdt  dt  4-5cc 


from  which,   with  eq    (IV-3), 


fcmiDCAA] 


Ar(app)  dt  dt 


(iv-15) 


cccc^BtAA]  A,^)+5Ccrcm)3:A.A:  AKilld)=  (DC)0.nd))oh3  +  S[cfc>vi)Ji:AA]  (DC)M 


)long 


Further  substituting  eqs  (IV-1)  and  (IV-2)  into  eq  (IV-16), 


_J  [_  ^Cc  (cm  Tic  a  A] 


.  a  j  JcacmncAAi 

'"Vtrc 


roe) 


'[cfcknrjCAA] 


W2    5 

NE    wCc(cmuLA'A] 


A,  ftrj=<OAP+  ^'""'  (c)Af+  |^fl£«i(-C)A  + 


COAr 


[cCc-cvOKa'A] 


tc(cmUC*AA] 


'Cc(cm,)]L'aa] 


(IV-16) 


(lv-17) 


The   lowest   order  of  the  forced  error  in  this  control  equation 
is  that  of  the  acceleration  of  the  wind.   The   coefficient  of  this 
term,    however  is 


sfc(CM)I  EAA] 
LSCC(CM)][AA]        W2NE   S[c(CM)][,AA] 


.  If 


sfC(CM)]  [*A*A]  is  chosen  as  equal  to  -1 —  ,  the  acceleration  .term 

w2NE 

disappears  and  the  control  equation  indicates  that  the  lowest 

order  of  forced  error  is  that  resulting  from  the  time  rate  of 

change  of  the  acceleration.  The  coefficient  of  this  term  is  of 

similar  form  to  that  of  the  acceleration  term,  and  is 


*N 
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1- 


S[C(CM)1    [AAl 
W  NE  S[C(CM)]    [AA]. 


It  thus  appears  that  the   acceleration 


rate  term  can  be  made  to  disappear  if  LC \ CM u  L  AA       ia   cnosen 

b[C(CM)][  AA] 

as  w  NE  * 

From  Routh's  stability  criteria  for  a  cubic,  the  sign  of 

each  term  in  the  homogeneous  equation  must  be  Dositive,  and 


CcfcmVlCAAl 


5 


[ckmllLAA] 


[cfcmMAAl 


[c(cw)]['aa] 


> 


[CCCW)]LA  A] 


(lv-18) 


Equation  (IV-18)    shows  that,   with  S  [c(CM)][  *A*A]   "  ~~2 —   »    the  system 


W 


NE 


is  stable  only   if 


;[C(CM)]  CAA] 


is  greater  that  W^ 


NE 


The 


S[C(CM)][AA] 

acceleration  rate  terra  can  therefore  not  be  eliminated  if  the 
system  is  to  remain  stable.  The  equation  that  has  been  considered 
here,  if  it  is  physically  possible  to  mechanize  it,  appears  to 
provide  a  satisfactory  system.  It  remains  to  be  seen  whether 
additional  terms  in  the  mechanization  equation  can  improve  the 
response. 

E.   Consider  next,  therefore,  the  oerformance  which  results 
when  to  the  mechanization  equation  just  considered,  the  first 
integral  of  Ar(app)  is  added  to  the  second  and  third.  This  mech- 
anization can  be  written 


JJ/Km  ***+V*ob *j?Ki*  dt  * +sbwiow /*►«.*  dk-t-StUKaaaArlM)=JffK^  dtdtdt  +5tc,cm)]t,w^/A„ppl  dt  dt  +\m^,u,Jk^  *    (IV-19) 


CONFIDENTIAL 


from  which,  with  eq    (IV-3), 


Lc(cm)]ffAj 


brUnM=    ^Sind)|onq     '    \fc«HlAAD    ^M,'^)  long    '       ^[C(cmUC  A  A]   (Dw(in<l)  Ug 


(IV-20) 


Further  substituting  eqs  (IV-1)  and  (IV-2)  into  eq  (IV-20), 


I      5 


CC(cw13C'aa1 


WNE     5lc<cw03 


H,  c+_,, ,\       i 


'[CfOYlflLAA] 


(trot) 


[AAH 


W2      5 

NE  J[C(C^][AA]/ 


A  ,       + 


>)-(tru.e.) 


**ME         "^ccfcm^LA  a] 


K(^T^K+  |^^(t)A  r+^=^L(t)Ar+ 


c-  ^/'  V        c 

^[cfcwBCA  AT  ^ 


[c(cml]L"A  A] 


[c(cw)][a'A] 


(C)Ar 


(IV-21) 


Here  again,    the   lowest  order  forced-error  term  contains  the 
acceleration  of  the  wind.   The   coefficient  of   this  term  is   now 

,     and  does  not    contain  any  controllable   sensitivity.   This 


W 


NE 


system  is  therefore  less  satisfactory  than  that  of  par  2-D. 

F.     For  the  sake  of   completeness,   Ar(app)   is  next  added  to 
its  integrals  in  the  mechanization  equation  just  discussed.    This 
gives 


AWmd)dtdtdt  +Scc(c^[AA3//Ar6hd)  dt  dt  +  5{c( HAM 


i  /  A,.6nd)  dt  -  pfcfcwaft'A]  Arfm<i)' 


from  which,  with  eq   (IV-3), 

m  •  •  *  •  • 

^tind)lo«q~'"SfcrcmW[AA]  C  Wri^)l0^^[cfc^KAA3^^Cma)Wg"T  5tc(C*01CAAl  (DQmd)Umq  '  O 


ng   '    ^fc(cwO][AA]  ^  "     '(undjlo^  [cfcv 


Arapp)dtdtdt+5 


[c(cm11CAA3 


/ /Arcapp)  at  dt  +o[t,,cml]rXA]  /A  r^pp)  dt  -r  oCc(cm):i  C-A-  w  Ar(app) 
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(IV-22) 


(1V-23) 


CONFIDENTIAL 


Further  substituting  eqs  (IV-1)  and  (IV-2)  into  eq  (IV-23) 


w 


^Vftrue) 


I       5, 


NE 


W*   C         /"Vtru«)    \Ji    e         "vftrue)   \,/t 


'ME  &fcin)JLAA] 


I  5 
T  5 


I 


W1  5 


A.*....,  =(C)Ar+  S(mm'xu\C)Ar  +  ^^L(c)K  + 


-ft  rue) 


'rc(cm)][AW 


A] 


5 


fc)A, 


[cfcn-OH  AA] 


(IV-24) 


This   system,    like  that  of   pa  r  2-E,    contains  an  unremovable 
steady-state-  error  proportional  to  the  magnitude  of  the  wind 
acceleration. 

G.     When  examining  possible  mechanization  equations,    it 
next  seems  reasonable  to  add  to  the  system  of  par  2-D  the  fourth 
integral  of   Ar(ind)»    us^ng  tne  fourth  integral  of  Ar(app).    This 
gives 


A^dtdtdtdt  +  S, 


rcfcmyjtAA.] 


ArMdtdtdt+5 


[cfcmyCA 


JjK<\rt<k&+\tcrtn*jKtoi)  *+V*)3tXiaArtI-ni,=  /WAr^p)dt  dt  dtdt 


from  which,  with   (IV-3), 


^[C(cm)][AAl  ^r(ind)   '    ^CcfcmutAAl  "r(ind)  +  ^fc(cm\][A>]^Hind)   '     ^[cfcmul'A"A]  'Mi'nd)- '  ^     'Cindll 


Further  substituting  eqs    (IV-1)   and   (IV-2)    into   (IV-26), 


(IV-25) 


(IV-26) 
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Arftr 


5        is 

^V  ft  rue) 


CcCcitJKaa] 


5 


•^OCa'a] 


CcCcm^rAA] 


WIS 


NE    ^tCfcW 


— k^-g""0^  A^=(OA**+fas-a'1  <£)A,+^^  fC)Ar  + 


Cc(cm)]  Caa3 


(C)Ar  + 


Cc(cwnCA~A3 


*u 


fOAr 


(IV-27) 


(ci-tOK'a'aJ 


As  was  true  with  the  corresponding  system  of  par  2-C,  an 
error  caused  by  the  velocity  of  the  wind  results. 

H.  When  the  third  integral  of  Ar(app)  is  added  to  the 
fourth  integral  used  in  par  2-G,  the  equation  becomes 


A^,  dt  dt  dt  dt  +  SK(cmaUA]  ///Ar(inJ)  dtdt  dt  +  S 


from  which,  with  eq   (IV-3) 


CCfcmtfCAA] 


Arffnj)  dtdt43DC(tmac.A.A1  //Vfinj)  dt+5cc(cm,]rXA3Ar(lhJ)=     I 7/Lpp)  dt  dt  d:  dt  +  S 


[cfcrrOJO 


JJjArh?p,<kdtdi 


^CdcmlKAAl   A(md)~'~^Lc(cw>)D[AA3  ^rCind)         •'WmncTA]  ^V(md)        ^^'fi'hd)lon4    '     ^ccCarOKAA}  VlAVfmd 


Further  substituting  eq    (IV-1)    into  eq    (IV-29), 


S  5 

rcfcwXlfA'AT  VVN£,-,Ccfcm1ilt'rAlJ 


"  fcrtl«) 


'CCfcm-PCAAl 


'cc^cmnfA-A] 


I 


W1   5 


r  A  ft  .= (C)Ar + 3cfcw,3cXA]  fc)Ar-h  -^aaa  fc)^+  ^^  ^  ^       I — (c)^ 


Acfc 


m)3C*AA3 


Ccfcmqc'AA] 


'ccfchOH'X'A] 


Ccfem!  KA  A3 


(IV-23) 


(IV-29) 


(iv-30) 


Here,   as  with  the  system  of  par  2-D,   the  control  equation 
indicates  that  the  lowest  order  of  error  in  the  steady  state 
results  from  the  acceleration  of  the  wind.   The  coefficient  of 
thiatenais     StC(CM)][AA]  1 _       . 

s[c(cm)][Xa]      v2m  s[c(cm)Haa! 
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cc 

As  before,    if  S  [c(CM)l    [*AA]    is  raade  equal  to        a        ,    the 

w  NE 

acceleration  term  disappears.    The  coefficient   of  the  acceleration 

rate   term  is        £MCM)UJAl  _  SLC(CM)]UA] .     #      „.   appears 

S[C(CM)]UA]       w2ne  s[c(cm)K7a] 

that   this  terra  can  be  made  to  disappear  if        [C  ( CM )]  [  AA 1        ^s 

s[c(cm)][aa3 

chosen  equal  to  J-   , 
W2NE 

From  Routh's  stability  criteria  for  a  quartic,  the  3ign  of 
each  terra  in  the  homogeneous  equation  must  be  positive,  and 


O  0  e 

^rc(cm):c*AAl    ^[C(CW)]tAA]     °[c(cm1][AA]  ^[cfcm)] 


LA  A]   pCcfcmKAAD    J  l^[cfcm13LA'*]J 


If  S[c(CM)][AAl  is  made  eclual  t0       ?   ■     i    eq    (IV-31)   becomes 

W  NE 


>0 


( lv-31 ) 


s 

!°[cfcm)]:/>  A] 


'rc(Cfn)3CAA] 


TcCcmDCAAJ 


w 


ME 


5  ^5 

-'CcfcmUQA  AD  J  ^[c  (cm)  3  [a  A] 


'[cfcm)][  a  A] 


5 


A  A] 


J[c(cmi]f 


A  A] 


>o 


(iv-32) 


Therefore,  in  order  to  eliminate  the  acceleration  rate  term, 

"  S[C(CD)][Ta1  [s[C(CU)1  Caa]]2  mu3t  be  grater  than  zero,  which 
is  possible  only  if  S  [c  ( CD  j]  [  *AA  ]  becomes  negative.  Routh's  first 
criterion  of  stability,  which  forbids  variation  in  the  signs  of  the 
quartic  coefficient,  is  then  violated,  and  the  system  is  unstable. 


rr 


ONFtDE;;, 


Once  again,  therefore,  it  has  not  been  possible  to  eliminate, 
in  the  steady  state,  any  effect  of  wind  of  higher  order  than 
acceleration.  Mechanization  eq  (IV-23)  appears  to  give  the 
same  accuracy,  in  the  steady  state,  that  results  from  the  use 
of  eq  (IV-15).  In  choosing  between  these  equations,  it  is 
necessary  to  compare  dynamic  response  and  ease  of  mechanization. 
Before  an  examination  is  made  of  this  problem,  a  study  will  be 
made  of  other  mechanization  equations. 

I.   In  an  effort  to  eliminate  the  effect  of  acceleration 
rate,  the  second  integral  of  Ar(app)  is  next  added  to  the  mech- 
anization equation  of  par  2-G.  This  gives 


JJ/fr,^  dt  dt  dt  dt  +  \a,^.JJfA  „,„,,  dt  dt  dt  +SM<^a,  fjh,  M1  dt  dt  +SWlwr.,y\„W)  A  dt  +3ow«,  A,,,.,  -fjJfK>„  *  dt  ck  dt  +SMnJffA,llp,  dt  dt  dt  +•  3^ 


CAA1 


from  which,   with  eq    (IV-3), 

Further  substituting  eqs    (IV-1)   and   (IV-2)    into  eq    (IV-34) 


)  long 


Km)  dt  dt 

(iv-33) 


(iv-34) 


—  )  |_  ^Cc(cm13CAA1 

1   wl  s 

'«£  VCC(cm)]c'AAl 


"r(true) 


ICCfcmOtAAj 


5 

^ccicmac'A'A:  *^N6  ^cc(c>tila'a: 


A,  (trut) 


w1  s 

NE  ^Cc(cm)][  A  A} 


Ar(truei  =£)A>  |^^  fc) Ar  +f^^  (c)Ar+  ^«^ 


^Cck^\]CAW 


'Cc(cm)3C'A"Xl 


5. 


(OAr 


:cic»m  a  a] 


I 


^tctcmniAA} 


fC)A, 


(iv-35) 


This  control  equation  shows  that,  instead  of  improving  the 
response,  the  use  of  eq  (IV-33)  has  resulted  in  an  acceleration 
error  which  it  is  no  longer  possible  to  eliminate. 
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J.   A  final  effort  to  remove  the  acceleration  rate  term  is 
made  by  adding  the  fifth  integral  of  Ar(app)  and  of  Ardnd)  to 
the  mechanization  equation  used  in  par  2-G.  This  give3, 


Arr,nd)  dt  dt  dt  dtdt-S^^Jjff^  dt  dt  dtdt+5cc^)n,Aw^AK,nd)  dt  <*t<it+ S^^ffA^dt  A+S^^A^A  +  Sukmmn  K^=[[[[[k^A  <*  dt  dt  dt  +S^lkJ[f[Aran)  dtdtdtdt 


from  which,    with   eq    (IV-3), 


^Cc(c^][XA3^i-<i»><i)  Ccfc~UtAAl^r(,«a^  '    *^Cc(cv,OtrA]  "KinaV^Wfc^tacXAl  "f(i»<IJ       '  ''Hind)  long    ""  ^  Cc  (c  mrj  c  AA3    lV\-)(ini)Uli 

Further  substituting  eqs    (IV-1)   and    (IV-2)    into  eq    (IV-37), 


(IV-37) 


-A, 


(W) 


wtcfc»«13LAA3 


(irut) 


5  5 

^LCtQvlH  A  AJ  ^CcfCmHCAA] 


*■*  CcfcmOfaVo  W  nE   ^Cc(cmUC"A"A]J 


A        —  < 


(tru«l 


,Cc(cmOLAAT 


'ccCcviUTa] 


wl  s 

WE         CCfcmnC   A 


"]  Arftrttf)=  fc)  A,  +  ******  (C)At+**™'Sn(C)K  +  ^^  (CK-f  ^^  CC)Ar  + ' (QAr 

'«)/  ^DifcmacTA]  ^Cc(cn,nc"A"A]  ^Ccfe»Tlt'X'«u  ^OtftmaCTw  *^[c(ewOn"X"A] 

llv-38) 


The  coefficient  of  the  wind  acceleration  term  of  the  forced 


«}■ 


error  now  has  the  form    K(CK)][  AAJ   _     1 

LsIc(cm)][Ta]         w  ne  s[c(cm)][Ta; 

To  eliminate  this,  SrcfcMlir  ^"1   must  be  made  equal  to  — * —  ,  as 

W  NE 

before.    In  this   equation,    the  wind  acceleration  term  has  the 


coefficient 


StC(CM)][AA]       _     3[C(CM)][AAl 

.  s[C(CN)][Ta]  W2NE  S[c(CMgf7A] 


If  this   is 


to  be  eliminated, 


EC  (CM)] 


_S[C(CM)1 


- I      must  be  made  equal  to     g» 

[AA]   J  W  NE 
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Using  Routh's  stability  criteria  for  a  quintic,  the  co- 
efficients of  the  homogeneous  equation  must  all  be  positive, 
and 


5      .  — 

K(cm)]CAA] 


tolltTAl 


tc(cmlH"AA] 


I      CcfcmmAA-]  CCfc^nCA'Al  \c(cWCAA)     ^CC(C^CAA]j 


WlKT«   ^^^{V*,,^     ^c^CAA3        ^mXAD  Acta*, C A A^ - {\<c»m? «]    {S, 


The  system  is  stable,  with  proper  choice  of  the  other  co- 
efficients,^ S  [C(CM)]L  A  A]   is  made  e<?ual  t0  ~r-   »  eliminating 

w  NE 

the  acceleration  term.    If  this   substitution  is  made  in  eq   (IV-39) 


and  if  in  addition  the  substitution        tC(CM)][  AA]     m         1 


tP(CM)][  AA] 


W 


is 


NE 


made, 


c«fcwatT4j^cefc«,)]rAK]       ic(c*m"£'nj 


>o 


(IV-39) 


°Cc(cm)][AAT  C 


[cfcm)3[A  A] 


CcrcmQu'A'Al 


w 


ME 


CcfcrnflC  X"  A] 


[CCc^DllA  A] 


ccfc^nc  AA] 


w 


N£ 


,yNE 


c>v)3C'a"a] 


V 


CcCc^DC'a'a] 


ME 


>0 


s  s 

^CcfcW)]£TAj   °tcfcw)Ji:AA] 


cc(c^):["aa] 


( IV-40) 


This  equation  simplifies  to 


80 


;:'i. 


CcfcmOCAA] 


CcCcm^c'A  a] 


[cfcnnOC  a'  a! 


w 


?       /  Wcc(c**03[a"a]   ^ 


.    -S 

CcCcm)3[AA]  Cc  Cc  mil  t  A  a] 


WE 


>0 


.w 


WE 


5  -  S 

^IcCc^llL'A  All      WCcfcm)3CA 


Al 


(IV-41) 


Since  all  of  the  sensitivities  are  positive,  this  equation 

2 


can  be  divided  by 


(wznrS[C(CM)][TAj' 

S[c(cm)][7a] 


,  provided  that 


only  S  [q  ( CD  )] ["aA  1     *3  maintained  smaller  than  (  A   ] 

\w  NE  / 

without  changing  the  inequality  sign.  This  yields 


5       ...   S 


[CfcmVJLAAj 


'[cternVJCVA] 

J 


Ml 


5, 


CcCcwl]  C  A  A3 


m 


CC(cnnlK'A"A 


]      -JccCc^iICaa]  ^[clfrnfl r  a" A]  J         J^   U 


( IV-42) 


or,   simplifying, 


-S 


£c(cm)l  C  A  A3 


>0 


Since  Sf;  ( CM )]  [*Ta  1   must  be  kept  positive  in  order  to  satisfy 
the  first  Routh  criterion,  both  criteria  cannot  simultaneously 
be  satisfied,  and  the  system  becomes  unstable  if  an  attempt  is 
made  to  eliminate  the  acceleration  rate  steady-state  error. 


(IV-43: 
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3 .  Implementation  of  Equations 

A.  Four  mechanization  equations  have  now  been  discovered 
that  eliminate   steady-state  error  caused  by  the  acceleration 
of  the  wind.  It  has  not  been  possible  to  eliminate  any  higher 
order  terms,  but  is  seems  unlikely  that  any  steady-state  error 
resulting  from  the  acceleration  rate  of  the  wind  can  have  appre- 
ciable effect.  The  simplest  of  these  equations  was  eliminated 
because  natural  frequency  could  not  be  controlled.  The  three 
usable  equations  are  (IV-15),  (IV-28),  and  (IV-36).  In  eq 
(IV-1S),  two  of  the  three  sensitivities  can  be  controlled;  in 

eq  (IV-23),  three  of  four  sensitivities  can  be  controlled;  and 
in  eq  (IV-3t>),  four  of  five  sensitivities  can  be  controlled.  In 
order  to  obtain  the  most  satisfactory  possible  dynamic  response, 
the  use  of  one  of  the  more  complicated  eauations  may  be  indicated, 
but  this  choice  of  equations  is  subordinate  to  the  considerations 
of  ease  and  accuracy  of  mechanization.  A  practicable  method  of 
implementing  eqs  (IV-15)  and  (IV-2S)  is  now  presented,  and  the 
complexities  compared. 

B.  In  mechanizing  eq  (IV-15),  (DC)  ( in(j)  (long)  is  u3ed 
instead  of  Ar(app)  as  the  output  of  the  pendulous  accelerometer, 
as  explained  in  par  1  of  this  chapter.  The  resulting  eq  (IV-16) 
is  restated  here  for  convenience,  in  integrated  form: 


A.,,.,,  dt  +  S, 


Art.M^/yhxtau.j  dt dt dt  +stet„w  JdcUh dt dt 


'ttCcm'OCA  a]    /  rV(m4)  u  ^tcCci^BtA  a] 

(IV-44) 


'J  8 


Now,  an  imaginary  auxiliary  direction  is  established,  as 
shown  in  Fig  IV-1.  This  direction  is  controlled,  with  respect 
to  the  reference  vertical,  by  the  equation 

A(auofreO—  ^EcfiioOUtAA]  /  /  (  W(ind)  \*n$  &   "^  (  |  V-45) 

The  controlled  member  is  positioned  from  the  reference  vertical 
through  the  use  of  the  auxiliary  direction  by  the  equation 


"Kind)      ^cc(c**)fau<y)[AA] 


//(Dc)(1-M<nuj  dt  dt-5Zc(CM)(i^[ilKJ  Ar(^Ahiikt  dt  ( iv-46) 


These  equations,  when  combined  with  the  relation 


(IV-4-7) 


become 


«/  </</  «///  uv-48) 

This  is  identical  with  the  desired  eq  (IV-44),  with  the 
relationships  among  sensitivities  that 


o  Jfcfc»)][AA] 

^[cfc^)fatm):[AA'J— 


^CCfc^lA  A] 


(IV-49) 
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ind)  ( long  ) 


"  (end)  (quk) 


Wcrt)fapp) 


(VertKn 


Cind) 


(Dir) 


f(in<4) 


(aux)  (ref  ) 


(Vert),ref 


(ref) 


(ayK) 


PIG.1Y-1 

AUXILIARY  DIRECTION  FOR  CUBIC  MECHf\NIIAT\ON 

EQUATION 
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Cc(ci^)(dUX)XAAl 


^[cfovO^CAA} 


IV-50) 


°[c(a 


[c  fa  i*y)"]  Ta  a] 


^[cCcwODL  a  AT 


Equations  (IV-45)  and  (IV-46)  can  be  mechanized  easily,  as 
shown  in  Fig  IV-2.  The  mechanization  of  eq  (IV-15)  is  therefore 
feasible. 

C.   Equation  (IV-23)  is  now  to  be  mechanized  in  a  manner 
similar  to  that  used  for  eq  (IV-15).  The  integrated  form  of  the 
equation,  using  (DC ) (ind ) (long)  instead  of  Ar(app),  as  before 
becones 


5        ...   A,-„+5 


[c  (c  **)]  c  a  a3 


Ahfl-^dt  +Sz«„mx>J  Arc^dt  dt=S 


lc(C^y]LAA]. 


(t)0(mJ)lo„,  dt  dt  dt 


For  this  equation,  a  pair  of  auxiliary  directions  is  esta- 
blished, as  shown  in  Fig  IV-3.  These  directions  are  controlled, 
with  respect  to  the  reference  vertical,  by  the  equations 


(IV-5D 


(DC)(iM)U^dt  dt  dt  dt 

(IV-52) 


Arau>vr^&~"5t:cfeux)aD[AAj//  (DC;(ihd)lo   dt  dt 
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(DC) 


ld)(OUX)Q 


(Lnd)donq) 


A(lnd)(Qux)b 


(Verthapp) 


(QuxKreDa 


(aux)(ref)b 


(Vert) 


(ref) 


(Vert)(Lnd)  {  ^r'(aux)a 

(Dlr>(QUX)b 


FIG.  IV-3.    AUXILIARY    DIRECTIONS    FOR  QUARTIC 
OR   QUINTIC    MECHANIZATION  EQUATION. 
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The  controlled  member  is   positioned   from  the   reference  vertical 
through  the  use  of  the  auxiliary  directions  by  the   equation 

^Mind)-         ^CcfGM&ioOHAA]    /™(i*d)(***)&  ^  ^Ccfc^faux):  Caa]// A(md)ra«.Ob    ^t  dt 

These   equations,   when  combined  with   the   relations 


(IV-54) 


(iv-55) 


(iv-56) 


(iv-57) 


become 


Sccfc^fewOHAAJ    l\(mi)  ut   "^"^ccfcmVauOKAAi  //AKl-Ma)  ^   "*     =SBAluft^CA'fa^CcCe»Vku<)3LAX]      //A^Qmd)!.^   ut    dt   Clt    ~H  O[c(aU)0b[AA]   O  ^^^[aa]  /  /  //  l^MmdJItmg    ^    "*    Clt    Ot 


Mt-(iVd)   ^Ccfe^lfau^DCA  AD  /  'VdVid) 

This  is  identical  with  the  desired  eq  (IV-52)  with  the 
relationships  among  the  sensitivities  that 


^CcCiux)a3LAA]         t 


'  Wcmu  C  A  A3 


(iv-59) 
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( IV-60) 


5, 


[c/crnVWODCAAj" 


CcfcmnCAAT 


CC«"i)]L'AA] 


(IV-61) 


'KtcmKauxnCAAJ 


CCfcrnQCAA] 


[c(cwO]tAAl 


(IV-62) 


Equations  (IV-53),  (IV- 54),  and  (IV-55)  can  be  mechanized 
almost  as  easily  as  were  eqs  (IV-45)  and  (IV-46),  as  shown  in 
Fig  IV-4.  The  choice  between  the  cubic  and  quartic  mechanization 
equations  is  largely,  therefore,  to  be  determined  by  a  comparison 
of  the  dynamic  responses  of  the  system  when  these  equations  are 
used. 

D.  The  method  of  mechanizing  the  quintic  equation,  (IV-36), 
is  similar  to  that  for  the  cubic  and  quartic  equations  just  dis- 
cussed. In  integrated  form,  using  (DC) (ind) (iong)  instead  of 
*r(app),  this  equation  becomes, 


W^-A-Aa  Ar^d)-1-  ^ccfcmtfc-A-Ai  /Ar6„d)  dt  +  Sccfc m):it.A. w  //Ar(lnd)  dt  dt ^5 c c (c w a cX A0  /  /  /Ar(md)  di  dt  dt  —  5 


Wcw)3Caa1 


foc)md),    dt  dt  dt  dt+//  /7(VC%^  dt  dt  dt  dt  dt 


(iv-63) 
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The  two  auxiliary   directions   established   by   eqs    (IV-53,    54,    56, and 
57)  are   required   to    implement   this   equation. 

The  controlled  member   is    positioned  by  the  use   of  these 
equations  and   the   equation 


When  these  are   combined,    the    control  equation  becomes 


(IV-64) 


Ar<;^H-Scc(cwX4U,)3l^% /Ar(md)dt  1  SCc(cm^au»mAAi  //Aw;nd)dt  dt  +JCcfcmKau*,3CAV] yy /Ar(md)dt  dt  dt      ^pei&u^b]tA»]  ^[ctcmHaurtitA-^ / / / /(l-'^Ld)io-1jdt  dt  dt  dt +oCc(aL)>)aKAA:s  o[cfeW)ftU<)][A-A]/////(DQ.nd) ,  dtdtdtdtdt 


(iv-65) 


This  is  identical  with  the  desired  eq  (IV-63).  The  mech- 
anization of  the  equations  from  which  (IV-65)  is  derived  is 
shown  in  Fig  IV-5.  This  mechanization  is,  of  course,  more 
complex  than  the  systems  previously  discussed.  Although  the 
mechanization  is  entirely  feasible,  the  quintic  equation  will 
not  be  considered  further  unless  the  dynamic  responses  of  the 
system  using  the  cubic  and  quartic  control  equations  show 
themselves  to  be  unsatisfactory. 
I».  The  Closed-loop  System 

In  order  to  examine  the  two  mechanization  equations  between 
which  choice  is  still  to  be  made,  it  is  necessary  to  examine  the 
response  of  the  closed-loop  longitudinal  system  to  wind  forcing 
functions,  using  each  of  these  mechanization  equations.  Response 
equations  for  particular  mechanizations  have  been  derived.  The 
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general  closed-loop  equations  are  derived  in  this  section. 

The  general  mechanization  equation  can  be  written,  for 
this  study,  as 


(iv-66) 


The  linear  distance  of  the  missile  measured  along  the  track, 
from  the  initial  reference  point,  can  be  written  as 


r~~  //  MS-fciro-dt  ut  +VDE-M3rmitiA|)t  +  ^cv,^d})  ( iv-67) 


Also, 

X 


~Ar(trae)  (IV-68) 


These  relationships,  with  eqs  (IV-1,  2  and  3)  and  the  small- 
angle  pendulum  equation 

•  •  • 

^C4ru01on3—        ^  ""  ^  (,V"69) 

Vie  ^iR 


combine  to  form  the   closed-loop  diagram  of  Fig  IV-6. 
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The  initial  conditions  in  eq  (IV-67)  do  not  affect  the 

e 

response  of  (C)Ar  to  V  rj;  m    (air)]  »  and  can  ^e  ta^en  as  zero. 
Figure  IV-6  can  then  be  simplified  to  Fig  IV-7. 

In  accordance  with  servomechanism  theory,  the  overall 
performance  equation  of  the  range  indication  system  can  be 
written 

rDc)   JcK  =  frFU4  frBLdo  -fl*L) 

^r/rris)  vv  I  \(pc)  (IV-70) 


where  (PF)(pcnd)   -  -i-    and  (PF)(kin)   -  -  -i-   .  The 

6IR  P  RE 

mechanisation  equations  to  be  examined  can  be  written 


I  ~l~   5CcfcwOCAA3  P (IV-7D 


(PF) 

Cc<ci*03Ca 


(PF) 


"^Cc(c^)]LAA3  r 


rmccKKqturtW     c  hx-l-S  n*4-  S  ft4 

^Cc(cw,tlCAA3  r  ^fcfcmuftAl  P      '     ^Ccfc^)3CA'A3   T 

(IV-72) 

Using  these  equations,  numerical  values  can  be  assigned  to 
the  sensitivities,  eq  (IV-69) solved  for  various  types  of  wind 
forcing  functions,  and  the  responses  using  eqs  (IV-73)  and 
(IV-71)  compared.  Numerical  values  are  first  assigned  to  the 
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cubic  mechanization  eq  (IV -*1),  and  then  to  the  quartlc  eq  flV-72). 
5.   Selection  of  Numerical  Values 

A.    The  homogeneous  control  equation  for  the  cubic  system  can  be  fac- 
tored into  a  first  order  and  a  quadratic  term,  becoming 


1 


fr+zrOfp^CDIflWp+W^O, 


CT 


(IV -73) 


which  expands  Into 


2(DR)IV+ 


1 


(P) 


P*+ 


W2+ 


2(DR)W 

(ct;  J 


P+CT=0 


(IV -74) 


Comparing  coefficients  term  by  term,  with 

CCT) 


Cc(cm)JC'A'AT 


W 


(IV-75) 


CDR)  = 


W2-W 


N£ 


2W(CT)\A/:, 


(IV-76) 


Cc(cm)][AA] 


(CTJ+ 


2(PR) 
W 


(IV -77) 


Then,  for  any  desired  characteristic  time  and  quadratic  natural  fre- 
quency, a  damping  ratio  and  set  of  sensitivities  can  be  computed.    Notice 
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that,  once  the  characteristic  time  and  the  natural  frequency  have  been 
selected,  no  control  remains  over  the  damping. 

B.  The  homogeneous  control  equation  for  the  quartic  system  can  be  fac- 
tored into  two  quadratic  terms.  As  in  the  track  control  system,  these  are 
made  equal,  giving, 

(p2+2(DR)Wp+W2)(pl+2(DR)Wp+W>0  w^* 

This  expands  into 

p4+4(DR)W|D3+  2W2[2(DR)+l]pi+4(DR)Wp+\AT=0  ov-w) 

Equating  coefficients  with  those  of  the  quartic  control  equation,  with 

Sr^  /      virV"*A  i    =   — * —   in  order  to  remove  the  acceleration  term  in  the 
[C  (cm)]|A  A  J        w2 

w    NE 

forcing  function, 

S  = (IV-80) 

^Ec(cm)UAAl         va/4 


(DR)  = 


'W2-2W2 


Ne 


4W  (iv-*i) 

NE 


Srctcm)]r^1=4(DR)  av-«» 

W 
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Plots  follow  for  responses  of  the  cubic  and  quartic  control  equations, 
with  various  sensitivities. 
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CHAPTER  V 
CONCLUSIONS  AND  SUGGESTIONS 

1.   Disturbances  to  the  Control  Systems 

In  examining  the  success  of  the  mechanization  systems  whose  responses 
have  been  analyzed  in  the  preceding  chapters,  it  is  necessary  to  consider  the 
types  of  disturbing  inputs  which  will  affect  the  systems.    These  disturbances, 
or  winds,  are  of  three  major  types.    First,  there  are  gust  disturbances  of 
fairly  high  frequency  and  short  duration.    These  are  onormously   attenuated 
by  the  long  time  constants  of  the  systems  analyzed,  and  should  cause  no 
trouble.   Second,  there  are  acceleration  pulse  type  Inputs,  such  as  would  be 
encountered  when  crossing  weather  fronts,  where  the  wind  velocity  changes 
very  rapidly  to  a  finite  and  fairly  constant  value.   These  can  be  analyzed  by 
the  step  function  responses  of  Chapters  HI  and  IV.    In  Fig  V-l,  the  path 
through  the  "eye"  of  a  hurricane  has  been  approximated  with  step  functions 
of  acceleration,  and  in  Fig  V-2,  the  response  of  the  track  control  system  is 
shown  to  be  satisfactory,  even  with  this  violent  disturbance. 

There  is,  however,  a  third  type  of  disturbance;  that  in  which  a  wind  of 
approximately  constant  velocity  slowly  rotates.*    This  condition  has  been 
approximated  by  sinusoidal  Inputs  as  shown  in  Chapters  in  and  IV.   The  error 
if  the  sinusoidal  forcing  function  continues  for  a  long  period  of  time,  is  seen 
to  be  very  large.   A  wind  of  about  forty  knots,  rotating  with  a  period  of  about 
one  hour,  will  introduce,  in  the  steady  state,  deviations  from  the  track,  in 
the  track  control  system,  of  about  20  miles. 


*Ref.  A.  G.  Bogosian  Notebook  dated  March  18,  1947. 
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TRACK 


"EYE"  OF    THE 
HURRICANE 


VELOCITY 
—  I00MPH. 


100  MP  H. 


ACCELERATION 

2 


2^250  FT/MIN 


O 


22,500  FT/Mlty2 


jr_  _j^0,Z50  FT/Mirt 


FIG.  V-l.  APPROXIMATION   OF  A  HURRICANE   USING 

ACCELERATION    PULSES. 
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Three  ways  of  correcting  this  condition  suggest  themselves: 

a.  Build  systems  with  variable  coefficients,  which  provide  the  neces- 
sary high  damping  with  step-inputs,  and  have  a  different  response  to  such 
low -acceleration  inputs  as  sinusoidally  varying  winds. 

b.  Increase  the  order  of  the  track  control  system  so  that  the  forced 
error  results  only  from  high  order  greatly  attenuated  terms. 

c.  Modify  the  system  to  reduce  the  effects  of  the  forced  errors. 

These  three  methods  are  all  examined  in  the  next  section  of  this  chapter. 

2.   Suggestions  for  Improvement  of  System  Response 
A.     Track  Control 

a.  A  variable  coefficient  system  can  be  of  value  only  if  changing  the 
damping  and  period  of  the  system  will  improve  the  response.   Since  the  co- 
efficient of  the  acceleration  rate  forced  error  term  increases  rapidly  with 
increase  in  damping,  it  appears  that  this  should  be  possible.   A  series  of 
track  control  responses  using  the  same  magnitude  of  disturbance  but  dif- 
ferent damping  ratios  has  therefore  been  examined,  Fig.  V-3.    It  is  seen 
that  no  choice  of  damping  is  satisfactory.    If  damping  is  entirely  removed, 
the  forced  error  of  the  simplified  system  reduces  to  zero.   Actually,  a 
forced  error  will  still  exist  if  there  is  any  vertical  wind  acceleration,  or  if 
non-ideal  components  are  used.    Then,  if  the  frequency  of  the  rotating  wind 
is  equal  to  that  of  the  system,  very  large  errors  can  result.    This  approach, 
therefore,  does  not  provide  the  needed  improvement  of  response. 

b.  The  purpose  of  this  section  is  to  investigate  the  response  of  a 
quintic  performance  equation.    From  eq  m-26  the  performance  equation 
may  be  written: 
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PERIOD  OF  A(Vert)(tc)(max) 

(Tm)  ROTATING  WIND  ?N  MINUTES  OF  ARC 

(DR)  ROTATING  WIND  [EQUAL  TO  NAUTICAL 

MILES  OFF  TRACK] 

1.875  10  7.735 

1.875  15  14.310 

1.500  15  13.620 

1.050  15  7.600 

1.875  20  13.033 

1.875  30  9.548 

1.500  30  14.899 

1.050  30  13.652 

1.500  40  12.080 

1.050  40  14.096 


MAXIMUM  TRACK  ERROR  FOR  VARIOUS  (DR)s  WITH  ROTATING  WINDS 

OF  31.4  MPH 


FIGURE  V-3 
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For  this  study  it  is  assumed  that  the  denominator  of  the  performance 
equation  may  be  factored  into  the  form: 

(p  +  ^KpS  2(0R)W  p+  W*)2 


Let         Wo  KWuE  where       WNEA 


Re 


Then  expanding  and  equating  coefficients  of  like  items, 


»i//rxD\\A/                   1                     9lR _.  SlC(CD)][A  A)                                                                 v-2 

CT          ->LCCCD)][A  A]v[(aLv.)-M]  D[C(CO)][AA] 

2K    WNE+4K(DR)    WNE+— ^ W  NE  +  '5^^ 

2K-W-NE^4K2(DR)-W-NE+   4K3(DR)W-          =w.  5Kc)CCD)3tAAt       +     SracD)3CAA],                               v^4 

CT  S[(C)(C0)3C'A  "A!            J[C(CO)][AA] 

4K3(0R)Wne          K4W4         -   W\,  5CCCCO)][AA]                                                                       Y-5 

CT  3[C(CD)][AA] 
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K      WNE       .     W2      ^[CCCD)][A  A] 


CT  S[C(CD)][AA]  (V"6) 

Examination  of  eqs  (V-2,  3,  4,  5,  and  6)  reveals  that  if  K  and  (DR)  are 
selected,  then  CT  and  all  the  sensitivities  are  fixed.   General  expressions 
are  now  developed  for  determining  the  numerical  values  of  CT  &ad  ui  certain 
ratios  of  sensitivities. 

From  simultaneous  solution  of  eqs  (V-3  and  5) 

4K(DR)(K2-1) 


CT    = 


WNE(-K4+6K*-1)  (V"7) 


From  eq  (V-6) 


5CcCcd)J[aa]  K3W*E(-K4+6K*-1) 
5[cCco)]rAA]    ~  4(DR)(K2-1)  (v"8) 

From  eqs  (V_5  and  7), 

SLc(cd)][aa]  K2W^E(5K2-1) 

5[CCCD)][AA]     =  K2-l  CV-9) 


From  eqs  (V-4,  6  and  7), 

7 


wNE 


S[C(CE))][AA]       _ 
T~? - 

^[CCC0)3CA  A] 
K7+4  K  [3(DR)J-  2]+K3[l3+6(PR)2]-aK[l+2(PR)*] 


4(DR)(K2-1) 

uo 


(V-10) 
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From  eqs  (V~2,  4  and  6) 

9ir  _   |-K6  +  4KfeC2-3CDR)2J+2K4[-7+4(DR)e]~4K'C2-3(DR)23-l  1  V_li 


'CCtCD)]['A  AJ  V[(Qir)-M] 


4K(DR)(K2-1) 


Limited  time  permits  the  evaluation  of  only  one  numerical  solution  of 
the  performance  equation. 
Let  K  =  2  and(D*)=  1. 


Then, 

"[C(CD)][AA]  x^    »VNE 


5  rrrr  mi  ta/a  T  i4   l/v. 


S[C(CD)][A  A] 


S[C(.CD)HA  A]        _     76WNE 


5[C(CD)]CA'A] 


^[C(CD)DCAA]      __  103  WNE 


'[CtCD)][A  A] 


<3iR  625    w 


5[C(CD)][AA]V[(Qlr)-IVl]  24 

1      _      7 
CT  24 


NE 


w 


NE 


Substituting  these  values  into  the  performance  equation  and  using  as  an 
input  Vr£  _  /ai  vi   =   2000  ft/min2  sin  —  t, 
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(C)  A  (Vert)  (tc)    =    .8  minutes  of  arc  (.8  nautical  miles) 
Comparing  this  answer  to  those  obtained  in  the  plot  of  results  in  Chap- 
ter in  indicates  that  the  quintic  performance  equation  produces  a  better  re- 
sponse than  the  quartic  one.   It  should  be  mentioned  that  no  effort  has  been 
made  to  select  optimum  parameters  for  the  quintic. 

c.   It  is  seen,  by  an  examination  of  the  track  control  system  perform- 
ance function,  that  no  control  over  the  gain  setting  exists.    As  was  shown  in 
Chapter  in,  it  is  desirable  to  have  the  sensitivities  negative  for  the  quartic 
performance  equation.   Making  this  change,  eq  in-25  becomes; 
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w2 


(PF) 


(tc)(CL)[VA]    = 


TC(CD)1CAA]  VC(cur)-M] 
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'[CCCO^tAA] 


where  sensitivities  are  now  positive  numbers. 

If  a  controllable  gain  coefficient  could  be  introduced  into  the  numerator, 
so  that  an  attenuation  of  the  forced  error  could  be  provided,  then  the  system 
might  be  made  to  perform  without  excessive  errors.    This  could  be  done,  it 
appears,  if  the  wind  acceleration  term  could  be  attenuated  after  entering  the 
system.   This  cannot  be  done,  but  an  approximation  of  this  desired  result 
is  achieved  if  additional  feedback  loops  are  provided  as  indicated  in  figure 
V-4.    The  performance  equation  for  this  system,  when  (PF)(aero)(inH)  is 
assumed  to  be  exactly  equal  to  (PF)/ aero\|AA  1 
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(.PF)(meCh)=  MECHANIZATION  FROM  CHAPTER  THREE 


FIG.  V-4.  MODIFIED    TRACK   CONTROL    SYSTEM. 
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I  hen  the  derivative  auxiliary  feedback  sensitivity  is  set  at  zero,  this  becomes, 
ith  proportional  plus  first  and  second  integral  mechanization  equation, 


P-b \ s 1  + 

F^'itc)(CL)[vA- 


5[CCC0)][A^]     f5(QlAX)[ 


'LC(CD)KAA] 


i) 
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tCCCDlKAA]  9lR  \      ,        [SCC(COl][AA]  <-  M/Z       |  ^    ,    e  5CC(CD>HAA]   ...2  .-  S[c(C0,]QftA] _1    ...? 

e C ^ P    +1= —   +   3(aiu)CAA]  WNElp     +0(au*)LAA]    ■= WNE  p+0((W)LAA] WNE 

D[C(CD)J[AA]  ^rC(CDH[AA]  VC(q(.v-1-M]  /  W[c(CD|][.a  tfl  /  *>[C(CD)][AA]  ^LCICOUCAA] 

In  this  equation,  a  certain  amount  of  gain  control  has  been  introduced 
Into  the  numerator,  at  the  expense  of  introducing  a  term  giving  a  forced 
error  caused  by  wind  acceleration.    Routh's  stability  criteria  show  that  this 

system,  with  proper  selection  of  sensitivities,  is  stable.    Limitations  of  i 

time  did  not  permit  the  authors  to  make  a  complete  analysis  of  this  control 
system.    However,  a  sample  result,  selected  at  random  is  calculated.    The 
denominator  is  first  factored  into  two  unequal  quadratics,  and  the  coefficients 
squared.    They  become: 
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Plc(cd)][aa]    \    ,    «.  _vu«      ^  rw.2-j-w.2  +?mQ\.(riQ)_    WW,  V-16 


M-  14 


iScC(CD)3C 


— V  VotAAiWif  -[VV+VV+ZCDR^DR),    WtW2 

A  A]    / 


cc(co,]tAAi\  ^CDRhW1W^^-fCDR)2W1^W2 
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The  values  chosen  were 
S(aux)[AA]    =    0A 
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1.5W 
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2 
(DR)2 


'    WNE 
-    0.7 


These  gave,  using  the  equations  given  above,  (DR).    =    9.18, 
![CIC°"lflflll  =  0. 124 ,  (f  "="°>"»1  )  =    22.64  and 


IR 


rccco)]CAA] 


'[C(CD)][A(V] 


SrciCUKAA]^t(Oif)-Ml 


=  20.72 


The  magnitude  of  the  maximum  forced  error  was  found  to  be  28.5  nautical 
miles,  and  the  disturbance  which  gave  this  maximum  deviation  was  a  sinu- 
soidally  varying  wind  with  a  velocity  magnitude  of  265.5  nautical  miles  per 
hour.    This  is  a  considerable  improvement  over  the  uncompensated  system. 
Next,  S/       \f*  a  1  was  se*  a*  umtv  *°  remove  the  accleration  term,  and 


(aux)[AA] 


was  set,  at  a  random  value. 


The  performance  equation,  with  quartic  mechanization,  becomes 
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Equating  coefficients,  as  before,  to  those  of  a  quartic  with  the  desired 
characteristics,  and  selecting,  at  random,  the  following  values: 

S(aux)[AA]    =    20 

Wl    =   W2    =    2WNE 
Answers  were  obtained  making  r^ "1  ~  .0885 


TGCCOnCAA] 


(OR), *(0R)2  =  2.75 ,  ("CCCCP)]CAA]     =  478 andf      9lR 


\=3.25 


[CCCD)][AA]v[(ai»-)-M], 


With  the  same  wind  input  as  that  used  in  the  preceding  example,  the  maxi- 
mum error  was  found  to  be  48.8  miles.    This,  again,  is  an  improvement 
over  the  uncompensated  system. 

It  appears,  since  in  the  first  example  the  error  was  caused  chiefly  by 
the  acceleration  term,  and  in  the  second  example  by  the  second  time  rate 
of  change  of  acceleration,  that  a  combination  of  the  two  systems,  varying 
both  s(aux)[AAl  and  S/auxxr^Aj  might  give  a  satisfactory  solution.    Time 
prevents  such  a  study  here.    A  simulator  would  greatly  facilitate  the  work 
of  selecting  optimum  sensitivities. 
B.    Range  Indication 

Lack  of  time  prevented  an  elaborate  examination  of  the  response  of 
the  range  indication  system  to  large  ranges  of  inputs.   Since  the  track  con- 
trol system  indicates  that  the  critical  response  is  to  sinusoidal  inputs,  a 
series  of  solutions  was  made  for  cubic  and  quartic  responses  to  a  sinu- 
soidal acceleration  input  with  a  period  of  thirty  minutes. 

As  indicated  in  the  plots  of  Chapter  IV,  the  response  is  improved  as  the 
damping  ratio  of  the  quadratic  in  the  cubic  equation,  or  the  two  equal  quad- 
ratics in  the  quartic  equation,  is  decreased.    At  a  damping  ratio  near  0.3, 


116 


CONFi.       HAL 


3.  TRACK    CONTROL    MECHANISM     SUMMARY 


MECHANIZATION      EQUATION 


* 

AzCCm)      =    S[c(C0>]    [AA](DC)tl„dHtc) 


PERFORMANCE     EQUATION   AMD    REMARKS 


(C)A 


REw'NE 


(Vert)(tc)       <-  vy 

°[c(CD)][AA]V[ia,r)-M] 


^^(Vert)ac)  *  WNe(UACVert)(tc) 


w: 


ME 


<?  V  I  [E-  (air)] 

°[C(CD)][AA]    V[cair)-M] 


velocity   of  Wind  introduces  forced  error 


Af  (cm)  s  "  S[c(cD)][AA]y  fDC)c™»Htodt 


A         «-  S 

M2(Cm)         ° 


[C(CD)]rAA]/7(DC^od)(U)^  dt 


'^- '"(Vert)  (tc)  + 


w, 


ME. 


lL)A(^rt)(tc)   - 


w 


WE 


Vr, 


'IR 


'[c(CD)][AA] 


-n 


ir  ^°[c(CO)][AA]V[la>Y).M] 


THIS     SVSTEM    HAS    NO    DAMPING    AND     BECOMES      Ki  N  5  T  A  6  L  E     IF     ATTEMPT     IS     MADE 
TO    ELIMINATE      FORCED    ERROR      INTRODUCED     BY     A  CC  ELERPvTvON     OF  ThE     VMIMO. 


'^'"(Vert)(tO  "* 


°[C(CD)][AA]  V[<QiO-M]/     x« 


(C)A 


S  .V 

d[C(CD)][AA]T[fa>0-M] 


<C)A,v<rt 


'ie 


(VertXtc)  R_        T[E-(a;^,] 


ME-(a'.rV 


THIS.      SYSTEM     IS     UNSTABLE    BECAUSE   THE     FIRST    ORDER    "TERM     vs    MlSS\KJ& 


HII 


CONFIDENTIAL 


TRACK     CONTROL   MECHANISM    SUMMARY  (.COWT) 
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4.   Range  Indication  Mechanization  Summary 

Examining  the  control  equation  obtained  In  Chapter  IV,  it  is  seen  that 
it  is  possible  to  eliminate  forcing  function  terms  of  orders  up  to  and  in- 
cluding wind  acceleration  terms.    In  a  stable  system,  higher  order  terms 
cannot  be  removed.   A  general  performance  function  including  all  mecha- 
nization equations  that  eliminate  acceleration  and  lower  order  terms  can 
be  written 

ine)  S(n-1)T    p  +  5(n)T 

WcM(AA,(ionS)=    p"+5(1)Tp"-'+-.+  S(n_Z)Tp*     Or-ra 

! 

If  more  terms  are  included  in  the  denominator,  velocity  and  position 
errors  result;  if  more  teTms  are  included  in  the  numerator,  control  is  lost 
over  low  order  terms,  while  higher  order  terms  are  eliminated.    If  terms 
are  omitted,  instability  results. 

A  block  diagram  showing  a  method  of  mechanizing  the  performance 
function  eq  (vi9)  follows, fig  (V-5). 

For  both  track  control  and  range  indication  systems,  the  response  for 
sinusoidal  forcing  functions  of  84.4  minute  period  appears  to  be  the  same 
regardless  of  the  choice  of  damping  or  order  of  the  equation  in  the  - 
mechanization  system. 
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APPENDIX  A 
DISCUSSION  OF  SIMPLIFYING  ASSUMPTIONS 

Several  simplifying  assumptions  were  made  in  reaching  the  mechani- 
zation equations  and  performance  equations  discussed  and  analyzed  in  this 
thesis.    These  will  now  be  examined. 

1.    Omission  of  the  Geodesic  Acceleration  Term 

The  geodesic  acceleration  arises  from  the  fact  that  navigation  is 
performed  on  the  surface  of  the  earth,  where  the  gravity  equipotential 
forms  the  geoid.    This  geoid  may  be  considered  to  be  an  ellipsoid  with 
roughness.    The  deflection  of  the  vertical  due  to  this  roughness  never  ex- 
ceeds seventy  seconds  of  arc,  and  (especially  over  land  masses)  is  usually 
far  smaller. 

While  the  missile  travels  over  this  ellipsoid,  computation  assumes 
that  the  surface  is  a  sphere.    The  shortest  distance  between  two  points  on 
this  sphere  is  a  great  circle.    The  gravity  verticals  fall  in  this  great  circle 
plane.   The  great  circle  plane  transfers  to  the  surface  of  the  ellipsoid  as 
a  plane  section,  but  the  gravity  verticals  no  longer  lie  in  this  plane,  but  in 
an  arc  intersecting  it  at  the  points  of  departure  and  destination.    The  mis- 
sile, traveling  along  this  curved  path,  experiences  a  horizontal  accelera- 
tion.    The  radius  of  curvature  of  this  path,  which  is  a  minimum  at  45  degrees 
north  or  south  latitude,  is  never  less  than  one  million  miles.    The  geodesic 
acceleration  is  therefore  of  negligible  magnitude  for  missiles  with  veloci- 
ties of  a  few  thousand  miles  an  hour. 
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Two  figures  (Fig  A-l  and  2)  prepared  by  Dr.  W.  Wrigley  of  the 
M.I.T.  Instrumentation  Laboratory,  illustrate  this  geodesic  acceleration. 

2.  Compensation  for  Coriolls  Acceleration 

The  acceleration  of  the  Coriolls  arises  from  the  fact  that  the  mis- 
sile is  moving  in  earth  space,  which  possesses  an  angular  velocity  with 
respect  to  inertial  space.    This  is  demonstrated  in  the  accompanying  fig- 
ure (Fig  A -3).    Compensation  is  independent  of  the  heading  of  the  missile, 
depending  only  on  a  knowledge  of  the  earth's  angular  velocity  with  respect 
to  inertial  space,  which  is  accurately  known,  on  latitude  and  on  ground 
speed.    An  indication  of  latitude  and  of  ground  speed  are  available  within 
the  missile,  the  accuracy  of  which  are  dependent  on  the  success  of  the  con- 
trol system.    Errors  resulting  from  inaccuracies  in  these  measurements 
will  be  negligible  in  any  practicable  system. 

A  certain  amount  of  cross -coupling  of  the  track  control  and  range 
systems  will  occur  through  the  Coriolls  Computer,  but  it  will  remain  small 
if  the  component  of  missile  velocity  perpendicular  to  the  track  is  not  large. 

3.  Small  Angle  Pendulum  Assumption 

For  simplicity  of  calculation,  a  linearization  assumption  has  been 
made  in  the  output  of  the  pendulum.    The  pendulous  accelerometer  actually 
solves  the  equation 

a[E  -  (cm)](hor)     _   ^  (DC) 

em 

The  control  equations  are  based  on  the  assumption  that 

a[E  -  (cm)](hor)      _    ^^v 
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The  validity  of  this  assumption  has  been  tested,  by  comparing  the  out- 
put, for  the  track  control  system,  to  a  pulse  of  wind  acceleration  of  2700 
feet/min  ,  lasting  for  one  minute,  using  the  transcendental  and  the  simpli- 
fied equation.    The  results  are  shown  in  the  accompanying  figure  (Fig  A -4). 

If  necessary,  a  tangent  pick-off  can  be  used  on  the  accelerometer  in- 
stead of  the  linear  pick-off  postulated  in  this  discussion. 

4.    Neglect  of  Vertical  Acceleration 

If  the  missile  does  not  maintain  absolutely  constant  altitude  (or  if 
the  magnitude  of  the  gravity  force  varies),  the  pendulum  angle  will  not  be 
a  function  of  horizontal  acceleration  alone.  This  effect,  together  with  the 
effect  of  an  aerodynamic  lag  and  a  pendulum  lag,  is  considered  in  the  fol- 
lowing treatment  of  the  range  indication  problem.  There  are  changes  in 
some  of  the  geometric  equations,  as  indicated: 
Aerodynamic  lag  term: 


QtE-tCmvKhor)  1 


(A-l) 


Vt(Qir>-E]  (CT"Wro)  P  +  l 

Accelerometer  lag  term: 

(DC)(tru^  _  1  (A-2) 

( DO(trueKcorr)  (CT)(pU)P+l 


Ar(«^,     =     aie-^)Kho.->  (A.3) 
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(DC) 
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CtrueKCOrr)  qjR  -   aLE.(cm):(Vert) 


(A-4) 


The  mechanization  chosen  was  the  cubic  mechanization  equation  discussed 
in  Chapter  IV. 

The  control  equation  for  the  range  indication  system,  using  the  con- 
cepts outlined  above,  becomes, 


(C)A 


(9lR~   Q-LE-CcmlKVert))  //NTX                  2        ("«-                    ..       _                                                          _  <5Lc(cm)]LAA1  (C  T)(pq)    QIR               SlcCcm)UAA1     Q  IE- tcm)l(Vert)  (CT)(pu)   ]                  ~                    ...             Stctcni)KAAl  9lR      ,     S(XCCm»  (AAUE-CcmiXVert) 
p  ^  U(pu)   T  L°lClcm)](AA)       y1H  T  u(.E-Ccm)KVert)       "  p  *    p  1    r    T  °  LcCOXH  AA)        p +    p 


(.Vert )  _ 


VlE-(air)} 


[p3+S[C(cm>KAA]    p2  +SLc(cml]tAAl    p  +  S[c(cm)UAA1HqiR-aiE-(cn»KVItrt)ll(CT)(pU)  p  "*  1 1  l(CT)(aero>  p  +  l] 


(a-5) 


If  the  sensitivities  are  chosen  so  that  the  cubic  in  the  denominator  be- 
comes the  product  of  a  quadratic  with  (DR)  =  0.715  and  a  period  of  forty 
minutes,  and  a  first  order  term  with  (CT)  =  14  min,  eq  (A-5)  becomes 


IC)A(. 


Vert) 


[E-(cmi](Vern    f       1       (  ,   ,  Q  [E-(cnQ3(VeTfl    (rr)i  Q  [E- (cm)](Vet^  p     (CT)(pa)  /  .       Q  [E-(cm)](Vert)\ 

-q^ +L4WE(^      9lR  HC1V>— q^         3P    4W,E    I1  qIR        j? 
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The  effect  of  vertical  acceleration  becomes  very  clear.    If  the  mis- 
sile momentarily  loses  lift  and  "drops",  the  denominator  goes  to  zero. 
The  numerator,  meanwhile,  acquires  a  term  which  gives  a  forced  error 
from  wind  acceleration.   This  error  will  change  sign  with  change  in  the 
direction  of  vertical  acceleration,  but  the  effect  of  the  denominator  will 
keep  the  opposing  effects  from  balancing.   The  damping  in  the  denominator 
cubic,  however,  will  return  the  missile  to  the  track  after  the  vertical  ac- 
celerations have  passed. 

Even  with  a  downward  acceleration  greater  than  the  effect  of 
gravity,  the  system  retains  stability. 

Reasonable  values  for  (CT)(aero\  and  ^T^(du)  w*^  vary  *rom  a 
fraction  of  a  second  to  a  few  seconds.  These  will,  therefore,  have  little 
effect  on  the  system. 

5.    Assumption  of  Perfect  Aerodynamic  Response 

Paragraph  4  indicates  that  neglect  of  aerodynamic  response  terms 
can  have  little  effect  on  the  performance  of  the  longitudinal  system.    In 
the  track  control  system,  however,  with  a  fourth  order  equation  or  higher, 
the  system  appears  stable  if  a  positive  pendulum  angle  causes  either  a 
left  or  a  right  rudder  deflection.   It  seems,  therefore,  that  even  a  small 
aerodynamic  lag  might  cause  the  system  to  become  unstable.    The  problem 
was  therefore  recalculated  for  one  choice  of  (DR),  using,  instead  of  perfect 
aerodynamic  response,  a  pair  of  equations.    The  first  gives  the  desired 
heading  of  the  missile  as  a  function  of  pendulum  angle;  the  second  relates 
this  desired  heading  to  the  actual  heading  of  the  missile  through  an  aero- 
dynamic lag. 
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Az(cmxind)  ~  S[c(cd)][aa]  (DC)(in(n<tcr  S[c(co)][AA]y  (DC)(ind)(tc)(dt)  ~S[c(cd)][AA]  J  J  (DC)(ind)(tc)dt  dt  (A-?) 


A 


The  variation  of  the  response  of  this  system  from  the  idealized 
system  discussed  in  Chapters  IE  and  IV  could  not  be  detected  in  the  plots 
drawn  by  the  Differential  Analyzer. 

6.    Assumption  of  Perfect  Gyros  and  Accelerometers 

Since  this  system  is  an  angle  measuring  system  which  carries  as 
a  reference  an  indication  of  the  point  of  departure  dependent  on  the  accu- 
racy of  the  gyros,  any  drift  in  the  gyros  which  define  the  vertical  plane 
will  give  an  equal  angular  error  in  range  and  deflection.    The  third  gyro, 
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which  maintains  the  orientation  of  the  track  plane,  will  introduce  a  non- 
linear error.    The  gyros  have  been  assumed  perfect  primarily  for  sim- 
plification of  the  problem.   If  gyros  of  sufficient  accuracy  are  found  to 
be  impracticable,  the  necessary  inertlal  reference  can  be  found  from  ob- 
servation of  the  fixed  stars. 

Permanent  errors  in  the  zero  indication  of  the  pendulums  also 
introduce  an  error  on  the  earth's  surface  equal  in  angle  to  the  pendulum 
error. 

Single  degree  of  freedom  gyros  and  pendulous  accelerometers  of 
types  which  may  prove  satisfactory  are  being  developed  in  the  Instrumen- 
tation Laboratory  at  M.I.T.   These  units  are  shown  diagrammatically  in 
Figs  (A-6)  and  (A-7). 

7.   Omission  of  Constants  of  Integration 

Wherever  possible,  throughout  this  thesis,  constants  of  integration 
have  been  assumed  to  be  zero.    For  example,  in  the  range  indication  prob- 
lem, the  missile  has  been  assumed  to  have  no  airspeed.   This  does  not 
change  the  time  response  or  the  size  of  the  range  or  track  error  resulting 
from  wind  disturbances,  but  will  have  a  large  effect  on  the  distance 
travelled  while  the  system  is  recovering  from  a  large  disturbance. 

At  all  times,  the  system  has  been  assumed  to  be  in  equilibrium 
when  hit  with  disturbing  wind  accelerations  in  the  time  solutions  of  Chap- 
ters in  and  IV. 
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APPENDIX  C 
GLOSSARY 

The  system  of  notation  adopted  for  use  in  this  thesis  has  been  selected 
because  it  fulfills  the  following  desirable  objectives: 

1.  It  is  easily  learned 

2.  It  is  adaptable  to  a  wide  range  of  situations 

3.  It  is  built  up  almost  exclusively  of  characters  found  on  the 
keyboard  of  a  standard  American  typewriter 

4.  Any  one  of  the  compound  symbols  of  the  system  is  readily 
interpreted  without  recourse  to  an  extensive  glossary 

A  few  simple  examples  suffice  to  explain  the  operation  of  the  notation 
system.    These  examples  are  shown  in  table  C~3.    The  short  table  of  key 
symbols  given  at  the  end  of  this  explanation  will  then  furnish  sufficient 
information  to  enable  any  compund  symbol  to  be  correctly  interpreted,  and 
to  provide- the  necessary  tools  for  the  generation  of  new  symbols. 

A  representative  list  of  the  primary  symbols  is  given  in  the  following 
table. 
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a  Acceleration 

A  Angle 

(CT)  Characteristic  time 

[(CT)PR]  Characteristic  time  -  period  ratio 

(C)  Correction 

(DR)  Dampting  ratio 

(Dir)  Direction 

(DC)  Dynamic  Correction 

(FR)  Frequency  ratio 

g  Gravity 

(PF)  Performance  function 

R  Radius 

S  Sensitivity 

(Sg)  Signal 

t  Time 

V  Velocity 

(Vert)  Vertical  direction 

W  Frequency 

PRIMARY  SYMBOLS 
TABLE  C-l 


* 
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(aero) 

Aerodynamic 

P 

Performance  operator  (d/dt) 

(air) 

Air  mass 

(pend) 

Pendulum 

(app) 

Apparent  value 

(pu) 

Pendulous  unit 

(aux) 

Auxiliary 

r 

Range 

(CD) 

Control  direction 

(ref) 

Reference 

(cm) 

Controlled  member 

(res) 

Resultant 

(dep) 

Departure 

t 

Time 

(dest) 

E 

(hor) 

(ind) 

IR 

Destination 
Earth 
Horizontal 
Indicated  value 
Inertial  reaction 

(tc) 

(true) 
(Vert) 
X 

Track  Control;  in  track  con- 
trol plane 

True  value 

Vertical 

Axis  through  nose  and  e.g. 
of  missile 

(kin) 
(long) 

M 

NE 

Kinematical 

Longitudinal  plane;  in 
range  direction 

Missile 

Earth  natural 

Y 
Z 

Axis  through  right  wing  and 
e.g.  of  missile 

Axis  through  e.g.  perpendicu 
lar  to  XY  plane,  directed 
downward 

MODIFYING  SYMBOLS  (SUBSCRIPTS) 


TABLE  C-2 
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Symbol 
V[E  -  (cm)] 

S[C(CD)][AA] 


(C)A 


(Vert)(tc) 


W 


NE 


AZ(cm) 

(DC)(ind)(tc) 

^(indjflong) 

Ar(ind) 


Definition 

Velocity  of  the  controlled  member  with 
respect  to  the  earth 

Sensitivity  for  correction  of  the  control 
direction  with  angle  input  and  rate  of 
change  of  angle  output 

Correction  to  the  angle  of  the  vertical 
in  the  track  control  plane 

Natural  earth  frequency  equal  to 

g 
IR 


*E 


Angle  about  the  Z  axis  of  the  controlled 
member 

Indicated  dynamic  correction  of  track 
control  system 

Indicated  dynamic  correction  of  longi- 
tudinal system  (i.e.  range  indication 
system) 

Indicated  range  angle 


Typical  Examples  Showing  How  Symbols  are  Compounded  from  the  Elemen- 
tary Forms 


TABLE  C-3 


141 


F  L     2  60 
A  P  I  3  6  i 


550-DUP 
BINDERY 

B  I  IIDERY 


^300G 


Thesis 

E$65         Erb 

A  theoretical  study  of 
automatic  inertial  navi- 
gation. 


FE     2  60     2 -*■ -C I        S5Q-DUP 
FEB   6   61  *SB3Tnv 


Thesis 

E565    Erb 

A  theoretical  study  of 
automatic  inertial  navi- 
gation. 
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